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the  broadband  method.  The  computer  programs  to  perform  these  two  analyses  are  written 
using  MATLAB  Examples  of  the  use  of  these  programs  are  included  in  this  report. 
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I.  INTRODUCTION 


In  recent  years,  the  rising  cost  of  machine  maintenance  has  driven  engineers  to 
develop  more  economical  and  efficient  methods  to  determine  machine  health  and  to  plan 
accordingly  the  required  preventive  and  corrective  maintenance.  The  most  p<^lar 
technique  in  use  today  is  a  predictive  maintenance  program  based  on  condition  monitoring. 
In  naval  applications,  condition  monitoring  is  commonly  achieved  utilizing  vibration 
measurement  and  analysis  on-board  surface  ships  and  submarines  [Refs.  1-7]. 

This  research  focuses  on  the  use  of  vibration  measurement  to  monitor  machine  health 
and  to  diagnose  system  problems  which  ccHild  lead  to  machine  failure.  In  addition  to 
providing  accurate  and  understandable  data  on  the  machine’s  current  condition,  a 
monitoring  and  diagnostic  system  must  also  limit  the  number  of  false  alarms.  Alarm 
threshold  settings  are  vitally  important  in  machine  vibration  diagnostics.  Alarm  thresholds 
set  too  high  may  result  in  premature  machine  failure  caused  by  an  undetected  failure 
condition  Alarm  thresholds  set  too  low  may  result  in  frequent  false  alarms  causing 
unnecessary  system  interruptions  and  repairs.  False  alarms  also  reduce  operator  confidence 
in  the  monitoring  and  diagnostic  system. 

Because  the  optimum  setting  of  alarm  thresholds  in  vibration  monitoring  and 
diagnostic  systems  continues  to  be  problematical,  the  goals  of  this  research  were; 

•  To  establish  a  statistical  analysis  method  for  setting  vibration  alarm  levels  by  using 
time  domain  data. 
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To  investigate  the  difference  in  broadband  alann  levels  between  the  cross  peak 


method  and  vibrometer  readings. 

•  To  perform  a  1/1  octave  band  analysis  for  setting  vibration  alarm  levels  by  using 
time  and  frequency  domain  data. 

Chapter  II  provides  the  background  information  on  the  machinery  utilized  in  this 
research  and  a  discussion  of  alarm  level  setting.  Chapter  III  describes  the  procedure  of  data 
acquisition  and  processing  and  the  methods  of  analysis.  In  the  time  domain,  a  statistical 
analysis  method  using  “cross  peak”  data  is  introduced.  In  the  frequency  domain,  a  1/1 
octave  band  concept  is  used  to  set  alarm  levels.  The  test  results  and  observations  using  time 
domain  data  are  presented  in  Chapter  IV.  Three  types  of  shipboard  fire  pumps  were 
analyzed  using  time  and  frequency  data.  Three  artificial  fault  simulations  were  performed 
to  assess  the  1/1  octav  band  method.  Chapter  V  demonstrates  the  1/1  octave  band  analysis 
results  using  the  NAVSSES  frequency  domain  database.  Chapter  VI  contains  conclusions 
and  recommendations  for  further  research. 


2 


II.  BACKGROUND 


A.  DATA  SOURCE 

The  fire  pump  on  a  naval  ship  provides  more  than  fire  fighting  water  to  damage  ccmtrol 
systems.  The  seawater  provided  by  fire  pumps  is  also  used  by  vital  air  ctmditioning  and 
chill  water  systems,  main  drainage  systems  and  decontamination  sprinkler  systems.  At  least 
one  fire  pump  is  running  whenever  a  ship  is  underway  or  at  anchor.  In  this  research,  twenty- 
three  fire  pumps  from  four  classes  of  ships  are  grouped  into  three  types  in  Table  1 . 


Table  1:  CLASSIFICATION  OF  FIRE  PUMP  DATA 


Hull 

Number 

Date  of 
Survey 

Number  of 
Fire  Pumps 

Number  of 
Pickup  Locations 

Number  of 
Data  Sets 

Type 

CG-59 

10/30/89 

6 

5 

30 

I 

CG-59 

2/26/90 

6 

5 

30 

DDG-994 

N/A 

3 

5 

15 

II 

FF-1062 

11/24/89 

4 

5 

20 

ni 

FF-1071 

3/10/89 

4 

5 

20 

The  vibration  velocity  signals  were  measured  by  means  of  transducers  strategically 
placed  at  pickup  locations  on  the  pumps.  Schematic  layouts  of  each  fire  pump  and  pickup 
locations  are  depicted  in  Figures  1  through  3.  A  list  of  abbreviation  for  pickup  locations  is 
shown  in  Table  2.  The  transducer  pickup  placement  is  uniaxial  with  one  radial  pickup  at 
each  bearing  and  one  axial  pickup  at  the  thrust  bearing.  This  arrangement  is  listed  in  the 
Vibration  Test  and  Analysis  Guide  (VTAG)  [Ref  8].  This  guide  provides  the  most  current 
technical  information  for  each  machine  in  a  particular  ship  class  and  includes  information 
on  pickup  locations,  operating  conditions  and  a  table  of  vibratitm  source  components. 
Table  3  through  Table  5  summarize  the  vibration  source  components.  These  tables  identify 
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exciting  elements  within  each  machine  and  list  the  vibration  frequencies  goierated  by  each 
element.  The  vibration  frequencies  were  normalized  as  multiples  of  the  machine’s  rotation 
rate,  called  “orders”.  The  time  waveforms  are  recorded  on  magnetic  tape  by  a  frequency 
modulated  (FM)  recorder.  Each  pickup  location  on  the  pumps  was  recorded  for  a  one 
minute  time  series  record.  Some  of  the  important  specifications  of  the  tape  rectx^der  are 
provided  in  Table  6. 

In  addition  to  the  data  tapes,  six  survey  dates  of  data  in  ASCII  form  were  received  from 
the  NAVSSES  database.  These  data  contain  the  frequency  spectrum  levels  measured  by 
vibrometer. 


TABLE  2:  ABBREVIATIONS  OF  PICKUP  LOCATION 


Abbreviation 

Description 

MB(FE) 

Motor  Bearing  (Free  End) 

MB(CE) 

Motor  Bearing  (Coupling  End) 

MB(CE/A) 

Motor  Bearing  (Coupling  End/Axial  direction) 

PB(FE) 

Pump  Bearing  (Free  End) 

PB(CE) 

Pump  Bearing  (Coupling  End) 

PB(CE/A) 

Pump  Bearing  (Coupling  End/ Axial  direction) 

UMB 

Upper  Motor  Bearing 

LMB 

Lower  Motor  Bearing 

LMB(A) 

Lower  Motor  Bearing  (Axial  direction) 

UPB 

Upper  Pump  Bearing 

LPB 

Lower  Pump  Bearing 
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Figure  1  Schematic  Layout  of  a  CG-59  Fire  Pump 


TABLE  3:  CG-59  VIBRATION  SOURCE  COMPONENTS 


Driven  (Pump) 


Driver  (Motor) 

Description 

Element 

Order 

Motor  Shaft  tRef.) 

1 

Fan  Blading 

5 

5 

Slots 

54 

54 

Bars 

44 

44 

Poles 

2 

2 

Bearing 

MRC  310 

Bearing 

MRC311 

Pump  Shaft 


Impeller  Vanes 


Bearing 


Element  Order 


1 


FAGWT 


MB(FE)  MB(CE)  MB(CE/A)  PB(CE)  PB(FE) 


Figure  2  Schematic  Layout  of  a  DDG-994  Fire  Pump 


TABLE  4:  DDG-994  VIBRATION  SOURCE  COMPONENTS 


Driver  (Motor) 

Driven  (Pump) 

Description 

Element 

Order 

Description 

Element 

Order 

Motor  Shaft  (Ref.) 

1 

Pump  Shaft 

1 

Fan  Blading 

7 

7 

Impeller  Vanes 

5 

5 

Slots 

48 

48 

E  :aring 

SKF306 

Bars 

38 

38 

Poles 

2 

2 

Bearing 

MRC310 

Bearing 

MRC311 
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Figure  3  Layout  of  a  FF-1062  /  FF-1071  Fire  Pump 


TABLE  5:  FF-1062  /FF-1071  VIBRATION  SOURCE  COMPONENTS 


Driver  (Motor) 

Driven  (Pump) 

Description 

Element 

Order 

Description 

Element 

Order 

Motor  Shaft  (Ref.) 

1 

Pump  Shaft 

1 

Poles 

2 

2 

Impeller  Vanes 

5 

5 

Bearing 

Ball  Bearing 

Bearing 

SKF6307 

TABLE  6:  SPECIFICATIONS  OF  FM  RECORDER 


Manufacture 

Dallas  Instruments  Inc. 

Model 

4800  FM  Recorder 

Frequency  Response 

2  to  SOOO  Hz  within  1  dB 

Signal  to  Noise  Ratio 

40  dB  rms 

Harmonic  Distortion 

No  harmonic  above  -40  dB 

Tape  Type 

TDK  AD  C60  cassette 

B.  Alarm  Level 

The  fundamental  question  that  must  be  answered  before  vibration  condition 
monitoring  can  be  used  as  a  diagnostic  tool  for  machinery  is  what  alarm  levels  identify 
operation-limiting  faults.  Generally,  the  vibrations  of  a  system  can  be  characterized  by  a 
reduced  data  set  in  various  domains.  The  criteria  used  to  set  alarm  levels  can  be  considered 
either  in  the  time  domain  or  the  frequency  domain,  as  discussed  below. 

1.  Time  Domain  Criteria 

a  Vibration  Severity  Criterion  Method 

The  simplest  time  domain  method  is  the  vibration  severity  criterion.  The  root 
mean  square  (RMS)  value  of  vibration  velocity  is  usually  measured  and  compared  with 
vibration  severity  charts  [Ref  9],  Figure  4  shows  the  vibration  severity  chart  that  was  first 
introduced  by  T.  C.  Rathbone  in  1939  and  then  later  refined  by  the  Instrument  Research  and 
Development  Corporation.  Various  national  standard  organizations  have  also  published 
standards  forjudging  vibration  severity.  For  example.  International  Standards  Organization 
(ISO)  standards  2372  and  3945  provide  severity  guidelines  for  machinery. 
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Figure  4  General  Machinery  Vibration  Severity  Chart  [Ref.  9] 


The  disadvantages  in  using  these  criteria  are  [Ref.  10]; 

•  The  criteria  can  be  used  only  for  specific  types  of  machinery  at  a  standard 

operating  condition. 

•  No  diagnostic  information  is  provided. 

•  The  criteria  are  less  sensitive  during  the  early  stages  of  damage. 

Thus,  the  vibration  severity  criterion  method  can  be  only  used  as  a  rough 
indicator  of  machinery  health. 
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b.  Amplitude  Probability  Criterion  Method 

A  More  sophisticated  method  was  developed  using  statistical  analysis  to  examine 
the  distribution  of  vibration  amplitudes.  The  amplitudes  used  can  be  either  peak-to-peak, 
peak  or  RMS  readings  of  displacement,  velocity  or  acceleration.  Both  Campel  [Ref  1 1  ]  and 
Murphy  [Ref  12]  use  this  statistical  method  to  establish  alarm  levels  based  on  the  mean  of 
the  reading  plus  3  standard  deviations.  The  key  drawback  of  this  approach  is  that  it  assumes 
that  a  Normal  distribution  of  the  linear  readings  exists.  Vibration  readings  have  a  “skewed” 
rather  than  “Normal”  distribution  [Ref  13]. 

2.  Frequency  Domain  Criteria 
a  Broadband  Criterion  Method 

The  broadband  criterion  method  utilizes  a  vibrometer  which  can  add  all  the 
energy  dissipated  over  a  wide  frequency  range  (typically  10  to  10,000  Hz).  The  overall 
energy  is  typically  calculated  by  applying  the  RMS  summation  method  to  the  spectrum.  If 
the  overall  energy  level  exceeds  a  predetermined  level,  then  alarms  are  triggered.  This  is  a 
simple  and  effective  method  in  most  cases.  But  in  some  cases  may  be  inadequate  because 
it  is  less  sensitive  to  small  changes  associated  with  bearing  defects  in  the  presence  of  large 
amplitude  frequency  components.  This  limitation  can  be  overcome  by  using  the  octave 
band  or  narrow  band  method. 

b.  Octave  Band  Criterion  Method 

The  octave  band  criterion  method  is  often  used  in  acoustics  to  determine  the 
energy  level  changes  due  to  noise  and  vibration.  This  method  utilizes  a  constant  percentage 
bandwidth  to  divide  the  frequency  range  of  interest  into  several  bands  which  provide  more 
detailed  information  than  the  broadband  presentation.  A  commonly  used  bandwidth  is  the 
one-third  octave  band.  Early  researchers  used  this  method  to  check  the  change  in  each  band 
level  and  determine  if  the  amplitudes  exceeded  normal  values  [Ref  14-15].  Chapter  ID 
discusses  this  method  in  detail. 
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c.  Narrowband  Criterion  Method 

Since  the  broadband  and  the  octave  band  criteria  methods  lack  detailed 
vibration  information,  a  narrowband  criterion  method  is  gaining  popularity.  The  bandwidth 
may  be  up  to  10%  of  the  center  frequency  range  (typically  10  to  1,000  Hz)  The  improved 
resolution  is  generally  up  to  400  or  800  lines  over  the  frequency  range  of  interest.  With 
knowledge  of  the  forcing  frequencies  of  rotating  components,  the  narrowband  data  is  very 
useful  for  diagnosing  specific  faults.  This  diagnosis  can  be  time  consuming;  it  is  also  often 
difficult  to  explain  the  source  of  some  of  the  peaks  in  the  ^ctra. 

Frequency  criteria  are  used  in  U  S.  Navy  surface  ship  to  establish  alarm  settings  for 
machinery  [Ref  3].  Vibration  data  for  each  individual  machine  and  location  are  processed 
to  yield  mean  and  standard  deviation  values  from  frequency  spectra  data  Generally, 
broadband  levels  are  used  as  a  screening  tools,  and  narrowband  levels  are  used  as 
diagnostic  tools 
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Ill.  METHODS  OF  ANALYSIS 


A.  DATA  ACQUISITION  AND  PROCESSING  SYSTEM 

A  block  diagram  of  the  data  acquisition  and  analysis  system  used  is  depicted  in  Figure 
S.  The  Model  4800  FM  recorder  plays  back  the  machinery  vibratirai  data  tapes  to  generate 
an  analog  signal.  By  using  a  connector,  tins  analog  signal  is  distributed  to  the  data 
acquisition  software  and  the  oscilloscope  simultaneously.  The  oscilloscope  contids  the 

quality  of  the  data  by  monitoring  the  signal  time  waveform.  The  EASYEST™  LX  software, 
developed  by  Keithley  Asyst,  was  used  as  a  data  acquisition  system.  This  PC-based  data 
acquisition  system  allowed  us  to  digitize  analc^  signals,  convert  the  units  from  voltage  to 
velocity  and  store  the  velocity  data  sequences  into  the  hard  disk  of  a  personal  computer. 
Finally,  MATLAB  software  was  used  to  retrieve  these  data  sequences  and  perform  data 
processing  MATLAB  is  a  very  powerful  software  tool  for  signal  processing.  All  programs 
used  in  this  research  were  written  using  MATLAB.  The  synopsis,  (m-line  help 
documentation  and  program  code  used  can  be  found  in  Appendix  A.  Examples  for  the  use 
of  these  programs  to  perform  the  analysis  are  included  in  Appendix  B. 

All  of  the  data  sets  were  sampled  at  a  sampling  frequency  of  10k  Hz.  The  sampling 
duration  for  each  data  set  was  36.684  seconds. 
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Figure  5  Block  Diagram  of  Data  Acquisition  And  Processing  System 
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B.  TIME  DOMAIN  -  CROSS  PEAK  ANALYSIS 


It  is  obvious  that  the  peak  envelt^  distribution  of  vibratimi  time  domain  signals  is 
relevant  to  the  vibration  severity  of  machinery  component.  Vibration  severity  criteria 
which  uses  measured  RMS  values  to  determine  the  severity  of  component  vibration 
without  considering  the  dispersion  of  the  signal  is  truly  a  rough  guess.  A  better  ^>proach  is 
to  represent  the  severity  of  component  vibration  in  terms  of  a  percentage  accq>tance  level. 
The  percentage  acceptance  level  provides  the  percentage  of  outcomes  >^ich  will  not 
exceed  this  level  threshold 

For  a  Normal  distribution,  the  acceptance  level  is  closely  related  to  the  mean  (p)  and 
standard  deviation  (o).  For  example,  p+1.96a  corresponds  to  a  97.5%  acceptance  level 
which  means  only  a  2  5%  probability  the  signal  amplitude  will  exceed  the  p+1.96o  value. 
For  the  sake  of  simplicity,  the  most  popular  assumption  to  assume  the  probability 
distribution  of  the  signal  was  a  Gaussian  (or  Normal)  distribution  However,  for  rotational 
machinery  vibration,  since  we  are  only  concern  with  the  magnitude  of  the  signal,  the 
probability  density  function  has  a  skewed  rather  than  a  Normal  distribution.  At  this  point, 
it  is  worthwhile  to  review  the  statistical  moments  and  central  moments  of  a  random  variable 

(Ref  16].  For  the  discrete  random  variables  with  probability  Pt(x),  the  mh  moment  E(j^) 
is  defined  as: 


Pnx,)  (1) 

km  1 

The  first  moment  is  particularly  important  in  many  applications  and  is  given  the  name  mean 
value  (n); 

=  Eix)  Pr(x,)  (2) 

I 

Of  greater  significance  are  the  central  moments  which  defined  as; 
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(3) 


Eiix-vi,y]  = 

km\ 

The  central  moment  for  n=  1  is  zero.  The  central  moment  for  n =2  is  a  very  important  which 
quantity  called  the  variance  of  the  random  variable. 


£t(»-M.)’)  =  X***-**-'’ 

The  standard  deviation,  which  corresponds  to  the  dispersion  of  the  randan  variable,  can  be 
(Stained  by  taking  the  square  roa  of  the  variance, i  .e. 


o  =  J  Variance  (5) 

As  noted  earlier,  the  mean  and  standard  deviation  are  often  used  to  compute  the  acceptance 
level 

The  higher  order  central  moments  (/i>2)  are  often  normalized  by  dividing  by  the  mh 
power  of  the  standard  deviation.  The  third  and  fourth  normalized  central  moments  are 
primarily  used  to  indicate  the  shape  of  the  probability  density  function.  They  are  defined 
by  the  following  equations; 


Skewness  = 


(6) 


Kurtosis  = 


(7) 


The  skewness  shows  information  about  the  position  of  the  peak  density  relative  to  the  mean 


value.  The  kurtosis  indicate  the  spread  in  the  distribution.  For  a  perfect  Nomal  distribution. 


the  skewness  is  zero  and  the  kurtosis  is  three. 


Since  we  are  interested  in  the  peak  envelope  distribution  of  the  vibration  signal,  we 
use  the  cross  peak  data  points  instead  of  the  overall  peak  data  points.  The  cross  peak 
technique  can  extract  the  dominant  component  from  a  complicated  time  waveform  data  set 
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by  selecting  the  maximum  peak  value  between  evoy  zero  crossing.  Figure  6  illustrates  the 
difference  between  overall  peak  and  cross  peak  points. 

As  we  will  see  in  Chapter  IV,  the  probability  distributions  of  cross  peak  points  have 
highly  positive  skewness  and  can  be  approximated  as  a  lognormal  distribution.  This  means 
that  the  distributions  will  become  normal  if  the  data  is  transformed  by  a  log  cm*  natural  log 
algorithm.  The  method  chosen  to  accomplish  this  transformation  was  to  convert  the  linear 
velocity  reading  to  velocity  dB  (VdB)  reading.  The  VdB  is  defined  as: 

Velocity (VdB)  =  20  log^  (8) 

w 

where  I'^/is  the  reference  level,  normally  0VdB=10'*  m/sec. 

Figure  7  shows  a  flow  chart  for  our  cross  peak  analysis  technique.  First,  the  sampled 
data  is  imported  into  a  MATLAB  worksheet.  The  DC  offset  is  then  removed  by  subtracting 
the  mean  value  of  the  sampled  data.  The  cross  peak  points  between  every  zero  crossing  is 
then  found  using  a  subroutine.  As  can  be  seen  from  Table  1,  the  dynamic  range  of  the  tape 
recorder  is  40  dB.  This  means  that  amplitude  smaller  than  1%  of  the  maximum  peak  may 
be  distorted.  For  this  reason,  we  set  a  1%  threshold  to  eliminate  these  distorted  points. 
Equation  8  was  then  used  to  transform  the  linear  velocity  scale  into  the  VdB  domain.  The 

reference  OVdB  was  10'*  m/sec.  After  we  computed  the  mean,  standard  deviation, 
skewness  and  kurtosis,  plots  of  the  probability  density  function  were  generated.  Finally,  the 
statistical  analysis  report  can  be  sent  to  laser  printer. 
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Figure  6  Comparison  of  Overall  Peak  and  Cross  Peak  Points  (a)  Overall  Peak  Points 

(b)  Cross  Peak  Points 
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Figure  7  Cross  Peak  Analysis  Flow  Chart 
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C.  FREQUENCY  DOMAIN  -  1/1  OCTAVE  BAND  ANALYSIS 


This  method  divides  the  frequency  spectrum  into  constant  percentage  bands  having 
the  same  ratio  of  bandwidth  to  center  frequency  [Ref  17],  Each  band  has  an  upper 
frequency  limit  if 2)  and  lower  frequency  limit  (/}).  The  center  frequency  (/^)  of  any  band  is 
defined  as 


f,  =  ,0  m 

The  ratio  of  center  frequencies  of  successive  proportional  bands  is  the  same  as  f^fi  for 
any  one  band.  i.e. 


^1—^1—  fl  —  2"^^ 

/.  fc  ~  \/.  " 


(10) 


One  third  of  an  octave  is  defined  as  /t=l/3.  The  American  National  Standards  Institute 
(ANSI)  preferred  center  frequencies  and  pass  bands  for  1/3  octave  bands  is  tabulated  at 
Table  7. 

The  method  chosen  to  examine  the  frequency  domain  alarm  levels  in  this  research  was 
to  use  1/1  octave  bands  (w=l).  For  the  10-5,000  Hz  frequency  range  of  our  system,  the 
frequency  spectrum  can  be  divided  into  9  bands  using  the  ANSI  preferred  center 
frequencies  The  center  frequencies  and  pass  bands  covering  the  frequency  range  10-5,000 
Hz  in  1/1  octave  bands  are  given  in  Table  8. 
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Table  7:  ANSI  PREFERRED  CENTER  FREQUENCIES  AND  PASS  BANDS  FOR 

1/3  OCTAVE  BAND 


1  Center 

Pass  Band  (Hz) 

Frequency  (Hz) 

Lower  Limit 

Upper  Limit 

12  5 

12 

14 

16 

14 

18 

20 

18 

22.4 

25 

22.4 

28 

31.5 

28 

35.5 

40 

355 

45 

50 

45 

56 

63 

56 

71 

80 

71 

90 

100 

90 

112 

125 

112 

140 

160 

140 

180 

200 

180 

224 

250 

224 

280 

315 

280 

355 

400 

355 

450 

500 

450 

560 

630 

560 

710 

800 

710 

890 

1000 

890 

1120 

1250 

1120 

1400 

1600 

1400 

1800 

2000 

1800 

2240 

2500 

2240 

2800 

3150 

2800 

3550 

4000 

3550 

4500 

5000 

4500 

5600 

6300 

5600 

7100 

8000 

7100 

9000 

10000 

9000 

11200 

Table  8:  ANSI  PREFERRED  CENTER  FREQUENCIES  AND  PASS  BANDS  FOR 

1/1  OCTAVE  BAND 


Figure  8  shows  the  flow  chart  for  our  1/1  octave  band  analysis  using  the  MATLAB 
program  developed.  This  technique  used  frequency  domain  data  to  compute  the  l/I  octave 
band  levels.  In  order  to  have  better  resolution  at  lower  frequency,  it  was  necessary  to  use 
two  spectrum  data  with  different  sampling  frequencies.  For  the  time  domain  data,  the  lower 
sampling  frequency  spanned  2  kHz  and  the  higher  sampling  frequency  spanned  10  kHz 
Twenty  blocks  of  data  is  sampled  at  each  sampling  frequency,  a  block  of  data  being  2048 
data  points.  After  sampling,  a  FFT  was  performed  for  each  block  with  a  Hanning  window 
to  obtain  the  smoothed  linear  velocity  spectrum.  The  averaging  VdB  spectrum  was  then 
produced  by  transforming  the  linear  spectrum  to  the  VdB  domain  and  taking  an  average 
over  the  20  data  blocks.  For  the  frequency  domain  data  received  from  the  NAVSSES 
database,  10  order  and  100  order  data  sets  were  used  as  the  lower  frequency  spectrum  and 
higher  frequency  spectrum,  respectively.  The  1st  through  6th  1/1  octave  band  levels  were 
then  computed  using  the  lower  frequency  sampling  spectrum  and  the  7th  through  9th  1/1 
octave  band  levels  obtained  using  the  higher  frequency  sampling  spectrum.  Finally,  the  1/ 
1  octave  band  analysis  was  done  afrer  combining  these  two  ranges  into  9  band  levels. 
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IV.  RESULTS  OF  USING  TIME  WAVEFORM  TAPES 


A.  CROSS  PEAK  ANALYSIS 

1.  Probability  Distribution  of  Cross  Peaks 

The  available  data  was  digitized  and  then  grouped  into  diree  types  after 
examination  of  the  spectrum  patterns.  Typical  time  series  and  frequency  spectnuns  for  each 
type  of  fire  pump  are  depicted  in  Figure  9  through  Figure  1 1 .  The  data  cdlecticm  sheets  are 
summarized  in  Table  9  through  Table  13  for  each  class  of  ship.  One  set  of  data  is  missing; 
the  actual  number  of  data  sets  analyzed  was  134.  The  operation  speed  of  the  shaft  was 
selected  as  a  reference  speed.  This  speed  was  used  to  normalize  the  frequency  scale  into 
“orders”.  The  “B.B.  Level”  is  the  broadband  level  measured  by  vibrometer  when  the  data 
was  recorded.  This  level  represents  the  total  energy  of  the  velocity  reading  for  each 
measurement.  The  full  scale  represent  the  calibration  factor.  The  “B.B.  Aim  Level”  is  the 
broadband  alarm  level  set  by  NAVSSES  using  the  vibrometer  reading  data.  Both  the  “B.B. 
Level”  and  “B  B  Aim  Level”  was  measured  for  a  frequency  range  from  10  Hz  to  10,000 
Hz.  The  remarks  column  describes  the  quality  of  the  data  set. 

In  order  to  understand  the  distributions  of  the  cross  peak  points,  probability 
histograms  for  each  set  of  data  were  plotted  using  linear  velocity  scale.  These  results  are 
provided  in  Appendix  C.  Inspection  of  these  figures  shows  the  distributions  to  have  an 
exponential  shape  for  the  most  part.  With  a  linear  velocity  scale,  this  type  of  distribution 
means  that  the  population  is  mainly  dominated  by  small  peaks  which  is  an  expected  result 
for  machinery  in  good  condition.  Those  distributions  having  a  large  population  at  higher 
amplitudes  are  regarded  as  representing  damaged  pumps.  Examination  of  these  damaged 
data  sets  shows  their  broadband  levels  to  be  higher  than  the  others  and  closer  to  the 
broadband  alarm  levels. 
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VdB(OVdB=  ie-8iTi/sec) 


Figure  9  Typical  Time  Series  and  Frequency  Spectrum  of  Type  I  Fire  Pumps  CG- 
59(10/30/89)  Fire  Pump  Measured  at  MB(F£) 
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VdB(OVdB  =  1e-8nVsec) 


Time  Domain 


Figure  10  Typical  Time  Series  and  Frequency  Spectrum  of  Type  II  Fire  Pumps  DDG- 

994  Fire  Pump  Measured  at  MB(F£) 
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VdB(0VdB=1e-8m/sec) 


Figure  11  Typical  Time  Series  and  Frequency  Spectrum  of  Type  III  Fire  Pumps  FF- 

1071  Fire  Pump  Measured  at  UMB 
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Table  9:  DATA  COLLECTION  SHEET  FOR  CG-59  (lQ/30/89) 
FIRE  PUMP  TYPE  I 


Machine 

Name 

Speed 

(RPM) 

Fire  Pump 
#1 

3595 

Fire  Pump 
#2 

3595 

Fire  Pump 
#3 

3595 

Fire  Pump 
#4 

3595 

Fire  Pump 
#5 

3595 

Fire  Pump 
#6 

3595 

Full 

Scale 

(VdB) 


B.B.  Aim 
Level 
(VdB) 


Remarks 


MB(CE) 

102 

120 

109 

PB{CE) 

110 

120 

118 

PB(FE) 

111 

120 

119 

PB(FE/A) 

114 

120 

121 

MB(FE) 

97 

100 

107 

MB(CE) 

101 

120 

109 

PB(CE) 

108 

120 

118 

PB(FE) 

106 

120 

119 

PB(FE/A) 

107 

120 

121 

MB(FE) 

97 

100 

107 

MB(CE) 

100 

120 

109 

PB(CE) 

111 

120 

118 

PB(FE) 

109 

120 

119 

PB(FE/A) 

113 

120 

121 

MB(FE) 

100 

120 

107 

MB(CE) 

100 

120 

109 

PB(CE) 

108 

120 

118 

PB(FE) 

110 

120 

119 

PB(FE/A) 

111 

120 

121 

MB(FE) 

100 

120 

107 

MB(CE) 

101 

120 

109 

PB(CE) 

110 

120 

118 

PB(FE) 

111 

120 

119 

PB(FE/A) 

114 

120 

121 

MB(FE) 

99 

120 

107 

MB(CE) 

100 

120 

109 

PB(CE) 

113 

120 

118 

PB(FE) 

111 

120 

119 

PB(FE/A) 

111 

120 

121 

Ibble  10:  DATA  COLLECTION  SHEET  FOR  CG-59  (2/26/90) 
FIRE  PUMP  TYPE  I 


Machine 

Name 

Speed 

(RPM) 

Pickup 

Location 

B.B. 

Level 

(VdB) 

FuU 

Scale 

(VdB) 

BB  Aim 
Level 
(VdB) 

Remarks 

Fire  Pump 
#1 

3595 

MB(FE) 

98 

100 

107 

MB(CE) 

97 

100 

109 

PB(CE) 

106 

120 

118 

PB(FE) 

110 

120 

119 

PB(FE/A) 

111 

120 

121 

Fire  Pump 
#2 

3595 

MB(FE) 

97 

120 

107 

MB(CE) 

96 

120 

109 

PB(CE) 

108 

120 

118 

PB(FE) 

107 

120 

119 

PB(FE/A) 

108 

120 

121 

Fire  Pump 

#3 

3595 

MB(FE) 

99 

120 

107 

MB(CE) 

100 

120 

109 

PB(CE) 

116 

120 

118 

Damaged 

PB(FE) 

115 

120 

119 

Damaged 

PB(FE/A) 

117 

120 

121 

Fire  Pump 
#4 

3595 

MB(FE) 

98 

120 

107 

MB(CE) 

97 

120 

m^mn 

PB(CE) 

108 

120 

118 

PB(FE) 

111 

120 

119 

PB(FE/A) 

112 

120 

121 

Fire  Pump 
#5 

3595 

MB(FE) 

99 

120 

107 

MB(CE) 

101 

120 

109 

PB(CE) 

109 

120 

118 

PB(FE) 

110 

120 

119 

PB(FE/A) 

113 

120 

121 

Fire  Pump 
#6 

3595 

MB(FE) 

99 

120 

107 

98 

120 

109 

PB(CE) 

115 

120 

118 

Damaged 

PB(FE) 

112 

120 

119 

PB(FE/A) 

113 

120 

121 

Table  11:  DATA  COLLECTION  SHEET  FOR 
FIRE  PUMP  TYPE  D 


Machine 

Name 

Speed 

(RPM) 

Pickup 

Location 

B.B. 

Level 

(VdB) 

Full 

Scale 

(VdB) 

BB  Aim 
Level 
(VdB) 

MB(FE) 

97 

120 

108 

Fire  Pump 

3555 

MB(CE) 

99 

120 

112 

#1 

MB(CE/A) 

no 

120 

112 

PB(CE) 

no 

120 

117 

PB(FE) 

109 

120 

119 

MB(FE) 

91 

100 

108 

Fire  Pump 

MB(CE) 

99 

120 

112 

#2 

3555 

MB(CE/A) 

107 

120 

112 

PB(CE) 

108 

120 

117 

PB(FE) 

109 

120 

119 

MB(FE) 

102 

120 

108 

Fire  Pump 

MB(CE) 

103 

120 

112 

#3 

3555 

MB(CE/A) 

112 

120 

112 

PB(CE) 

113 

120 

117 

PB(FE) 

116 

120 

119 

Remarks 


Bad  Data 
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Table  12:  DATA  COLLECTION  SHEET  FOR  FF-1062  (11724/89) 

FIRE  PUMP  TYPE  in 


Machine 

Name 


Fire  Pump 
#1 


Fire  Pump 
#2 


Fire  Pump 
#3 


Fire  Pump 
#4 


Speed 

(RPM) 

Pickup 

Location 

B.B. 

Level 

(VdB) 

FuU 

Scale 

(VdB) 

BB  Aim 
Level 
(VdB) 

UMB 

110 

120 

108 

LMB 

111 

120 

110 

3570 

LMB(A) 

103 

120 

108 

UPB 

117 

120 

120 

LPB 

112 

120 

121 

UMB 

97 

120 

108 

LMB 

98 

120 

no 

3570 

LMB(A) 

98 

120 

108 

UPB 

107 

120 

120 

LPB 

105 

120 

121 

UMB 

101 

120 

108 

LMB 

104 

120 

no 

3570 

LMB(A) 

101 

120 

108 

UPB 

109 

120 

120 

LPB 

111 

120 

121 

UMB 

98 

120 

108 

LMB 

103 

120 

no 

3570 

LMB(A) 

98 

1 

120 

108 

UPB 

114 

120 

120 

LPB 

113 

120 

121 
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Table  13:  DATA  COLLECTION  SHEET  FOR  FF-lOTl  (3/10/89) 
nRE  PUMP  TYPE  in 


Machine 

Name 

Speed 

(RPM) 

Pickup 

Location 

B.B. 

Level 

(VdB) 

Full 

Scale 

(VdB) 

BB  Aim 
Level 
(VdB) 

Remarks 

Fire  Pump 
#1 

3570 

UMB 

98 

120 

108 

LMB 

99 

120 

110 

LMB(A) 

97 

120 

108 

UPB 

110 

120 

120 

LPB 

106 

120 

121 

Fire  Pump 
#2 

3570 

UMB 

96 

120 

108 

LMB 

104 

120 

no 

Missing 

LMB(A) 

102 

120 

UPB 

117 

140 

120 

LPB 

114 

140 

121 

Fire  Pump 
#3 

3570 

UMB 

105 

120 

108 

LMB 

102 

120 

LMB(A) 

101 

120 

108 

UPB 

111 

120 

120 

LPB 

107 

1 

120 

121 

Fire  Pump 
#4 

3570 

UMB 

96 

100 

108 

LMB 

102  ! 

120 

no 

LMB(A) 

99 

120 

108 

UPB 

114 

120 

120 

LPB 

109 

120 

121 

31 


2.  Statistical  Analysis  Results 

As  mentioned  before,  *he  distributions  of  the  cross  peak  points  in  a  linear  velocity 
scale  have  an  exponential  shape  histogram.  This  shape  histogram  can  be  approximated  as 
a  lognormal  distribution.  Because  a  velocity  decibel  is  a  unit  of  measurement  equal  to  the 
logarithm  of  the  ratio  it  is  convenient  to  use  the  VdB  transformation  to  convert  an 
exponential  shaped  distribution  into  a  Normal  distribution.  Before  we  applied  this 
transform,  the  1%  threshold  was  used  to  remove  distorted  data.  The  data  sets  with  same 
type  of  fire  pump  and  same  pickup  location  were  then  added  together  to  form  a  combined 
data  set  to  perform  the  statistical  analysis  The  results  are  collected  in  Appendix  D.  Figure 
D.  1  through  Figure  D.  1 5  present  the  statistical  analysis  results  for  all  of  the  available  data. 
Figure  D.16  through  Figure  D.21  show  the  statistical  analysis  results  without  including  the 
damaged  data  sets.  Table  14  tabulates  these  results  in  the  linear  velocity  and  VdB  domains. 
A  comparison  of  alarm  thresholds  between  the  broadband  alarm  levels  (B.B.  Aim.  Level) 
and  the  computed  p+2a  level  is  shown  in  Table  IS.  It  should  be  mentioned  that  these  two 
alarm  thresholds  are  based  on  different  methods  using  data  from  different  domains.  The 
broadband  alarm  level,  as  noted  earlier,  is  obtained  by  using  vibrometer  readings  in  the 
frequency  domain.  The  computed  p+c  levels  represent  the  97.5%  acceptance  level  of  the 
cross  peak  envelopes  in  the  VdB  domain.  Several  observations  can  be  made: 

•  The  distribution  of  the  cross  peak  points  in  the  linear  velocity  domain  have  a 

highly  positive  skew.  This  can  be  treated  as  a  lognormal  distribution 

•  By  using  the  Vdb  transformation,  the  distributions  in  the  linear  velocity  domain 
can  be  transformed  to  Normal  or  near  Normal  distributions  Generally  si>eaking,  the  data 
sets  at  pickup  locations  for  the  motor  have  a  skewness  less  than  0. 1  and  the  data  sets  at 
pickup  locations  for  the  pump  ends  have  a  skewness  varying  between  0. 1  and  O.S. 

•  As  can  be  seen  from  Table  14,  removing  the  damaged  data  sets  does  not  cause 
the  probability  distributions  in  the  Vdb  domain  to  become  much  closer  to  Normal  as 
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measured  by  the  skewness  and  kurtosis.  For  the  MB(CE)  pickup  location  of  a  type  11  pump, 
not  considering  the  damaged  pump  data  sets  makes  the  skewness  even  more  severe  than 
when  the  damaged  pump  data  sets  are  considered  Since  there  were  only  3  type  D  pumps, 
this  result  is  not  surprising. 

•  In  general,  the  B.B.  Alarm  Levels  and  computed  p-f2a  levels  in  the  VdB 
domain  are  quite  close.  The  relative  errors  between  these  two  methods  for  Type  I  fire 
pumps  are  less  than  S%.  It  is  also  noted  from  Table  IS  that  the  computed  p-f2o  levels  in 
cases  with  small  skewness  is  closer  to  the  broadband  alarm  levels  than  in  cases  with  high 
skewness.  This  is  because  the  computed  p+2a  levels  for  the  Normal  distribution  is  very 
sensitive  to  skewness. 
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Table  14:  STATISTICAL  ANALYSIS  RESULTS  FOR  CROSS  PEAK  POINTS 


Type 

Pickup 

Position 

Linear  Velocity  Domain 

VdB  Domain 

(mm/s) 

O 

(mm/s) 

Skew. 

Kurt. 

li 

(VdB) 

<T 

(VdB) 

Skew. 

Kurt 

I 

MB(FE) 

0.495 

0.533 

1.850 

6.201 

89.36 

9.20 

-0.046 

2.378 

MB(CE) 

0.539 

0.608 

2.052 

7.195 

90.20 

9.12 

0.049 

2.444 

PB(CE) 

2.828 

2.680 

1.517 

5.638 

104.6 

9.78 

-0.444 

2.476 

PB(FE) 

2.717 

2.681 

1.253 

3.959 

1036 

10.45 

-0.335 

2.232 

PB(FE/A) 

3.699 

3.307 

1.086 

3.515 

106.9 

9.90 

-0.519 

2.472 

2.319 

1.988 

1.046 

3.304 

103.3 

9.25 

-0.510 

2.506 

PB(FE)‘ 

2.525 

2.472 

1.169 

3.409 

103.0 

10.30 

-0.317 

2.228 

II 

MB(FE) 

0.504 

0.604 

2.059 

6.976 

88.4 

10.50 

-0.139 

2.420 

MB(CE) 

0.727 

0.763 

3.243 

66.02 

93.1 

8.78 

-0.022 

2.320 

MB(CE/A) 

2.148 

2.493 

1.321 

3.604 

100.0 

11.49 

0.019 

1.926 

PB(CE) 

2.484 

2.344 

0.905 

2.675 

102.6 

1084 

-0.368 

2.038 

PB(FE) 

2.166 

2.739 

2.008 

6.446 

100.8 

10.29 

0.113 

2.366 

MB(CE)’ 

0.619 

0.811 

2.146 

6.762 

90.5 

9.19 

0.495 

2.566 

III 

UMB 

0.672 

0.926 

4.985 

36.52 

92.1 

8.69 

0.029 

2.953 

LMB 

0.923 

1.046 

3.350 

19.59 

95.3 

8.56 

-0.077 

2.697 

LMB(A) 

0.703 

0.648 

1.881 

7.538 

93.3 

8.46 

-0.349 

2.663 

UPB 

4.141 

4.442 

1.570 

5.446 

106.1 

11.84 

-0.408 

2.207 

LPB 

3.160 

3.246 

1.779 

6.775 

105.1 

10.20 

-0.373 

2.430 

0.597 

0.602 

2.488 

11.46 

91.8 

8.34 

-0.167 

2.595 

HQHII 

0.776 

0.736 

1.980 

7.430 

94.3 

8.15 

-0.181 

2.555 

3.603 

3.906 

1.711 

6.361 

105.1 

1160 

-0.367 

2.176 

*  :  The  statistical  ana'ysis  results  without  considering  the  damaged  data  sets. 
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Table  15:  COMPARISOM  OF  COMPUTED  n+2o  LEVEL  AND  BROADBAND 

ALARM  LEVEL 


TVpe 

Pickup 

Position 

Skewness 

4 

(VdB) 

o 

(VdB) 

ll+2o 

(VdB) 

B.B.Alarm 

Level 

(VdB) 

Relative 

Error 

MB(FE) 

-0.046 

89.36 

9.20 

107.8 

107 

0.71% 

MB(CE) 

0.049 

90.20 

9.12 

108.3 

109 

-0.68% 

PB(CE) 

-0.444 

104.6 

9.78 

124.2 

118 

4.95% 

I 

PB(FE) 

-0.335 

103.6 

10.45 

124.5 

119 

4.42% 

PB(FE/A) 

-0.519 

1069 

9.90 

126.7 

121 

4.50% 

-0.510 

103.3 

9.25 

121.8 

118 

3.11% 

-0.317 

103.0 

10.30 

123.6 

119 

3.72% 

MB(FE) 

-0.139 

88.4 

10.50 

109.4 

108 

1.27% 

MB(CE) 

-0.022 

93.1 

8.78 

110.6 

112 

-1.25% 

11 

MB(CE/A) 

0.019 

100.0 

11.49 

123.0 

112 

8.91% 

il 

PB(CE) 

-0.368 

102.6 

10.84 

124.3 

117 

5.86% 

PB(FE) 

0.113 

1008 

10.29 

121.4 

119 

1.96% 

0.495 

90.5 

9.19 

108.9 

118 

-8.37% 

UMB 

0.029 

92.1 

8.69 

1095 

108 

1.40% 

LMB 

-0.077 

95.3 

8.56 

112.4 

110 

2.11% 

LMB(A) 

-0.349 

93.3 

8.46 

110.2 

108 

2.04% 

Til 

UPB 

-0  408 

106.1 

11.84 

129.8 

120 

7.54% 

111 

LPB 

-0.373 

105.1 

10.20 

125.5 

121 

3.59% 

-0.167 

91.8 

8.34 

108.5 

108 

0.44% 

E&33lll^l 

-0.181 

94.3 

8.15 

110  6 

110 

0.52% 

nnaH 

-0.367 

105.1 

11.60 

128.3 

120 

6.47% 

* ;  The  statistical  analysis  results  without  considering  the  damaged  data  sets. 
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B.  1/1  OCTAVE  BAND  ANALYSIS 
1.  1/1  Octave  Band  Analysis  Results 

The  1/1  octave  band  analysis  was  performed  for  all  types  of  fire  pumps  >vith  the 
time  domain  data.  Band  1  through  band  6  were  computed  by  using  fs-2  kHz  to  provide 
better  resolution  at  lower  frequencies.  Band  7  to  band  9  were  computed  by  using  fj^lO  kHz. 
Appendix  E  contains  figures  for  all  of  the  1/1  octave  band  levels  for  each  set  of  data.  The 
10-5,000  Hz  broadband  levels  are  also  shown  on  these  plots.  It  was  computed  by  adding  all 
the  VdB  spectrum  using  a  RMS  algorithm.  The  means  (p)  and  standard  deviations  (a)  of 
each  octave  band  levels  for  the  same  type  fire  pumps  with  the  same  pickup  locations  was 
also  computed.  Figure  12  through  Figure  14  show  the  summarized  1/1  octave  band  levels 
for  each  type  of  fire  pump  without  considering  the  damaged  data  sets.  The  thin  line 
represents  the  mean  (p)  band  level  and  the  thick  line  represents  the  mean  plus  one  standard 
deviation  (|i+<j)  band  level.  Table  16  tabulates  these  band  levels. 

Generally  speaking,  the  dominate  levels  are  located  in  the  first  six  bands  (10  to 
710  Hz).  This  implies  that  the  energy  of  vibration  is  concentrated  at  lower  frequencies.  For 
fire  pumps  with  a  3555  RPM  operation  speed,  these  bands  are  approximated  up  to  12 
orders.  Upon  examination  of  the  third  band  through  the  sixth  band  for  each  data  set,  it  is 
obvious  that  the  band  levels  for  damaged  data  sets  are  higher  than  the  for  undamaged  pump 
data  sets. 
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VdB  (0VdB=1  e-8  m/s)  VdB  (0VdB=1  e-8  m/s) 


1/1  Octave  Band  Analysis 
CG-59  Fire  Pump 

Each  Data  Set  Use  20  Averaged  FFT  Data 


(a)  Pickup  Location  :  MB(FE) 


Frequency  (Hz) 

(b)  Pickup  Location  ;  MB(CE) 


Frequency  (Hz) 


Figure  12  The  Summarized  1/1  Octave  Band  Levels  for  Type  I  Fire  Pumps 
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VdB  (0VdB=1  e-8  m/s)  VdB  (0VdB=1  e-8  m/s)  VdB  (OVdBsle-S  m/s) 


2 

3 

10 

10 

10 

Frequency  (Hz) 

Figure  12  The  Summarized  1/1  Octave  Band  Levels  for  Type  I  Fire  Pumps  (Cont) 
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VdB  (0VdB=1  e-8  nVs)  VdB  (OVdBsl  e-8  tn/s) 


1/1  Octave  Band  Analysis 
DDG-993  Fire  Pump 

Each  Data  Set  Use  20  Averaged  FFT  Data 


(a)  Pickup  Location :  MB(FE) 


Frequency  (Hz) 


(b)  Pickup  Location :  MB(CE) 


Figure  13  The  Summarized  1/1  Octave  Band  Levels  for  Type  n  Fire  Pumps 
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VdB  (0VdB=1  e-8  iWs)  VdB  (0VdB=1  e-8  nrVs)  VdB  (0VdB=1e-8  nVs) 


(c)  Pickup  Location :  MB(CE/A) 


Figure  13  The  Summarized  1/1  Octave  Band  Levels  for  Type  II  Fire  Pumps  (Cont) 
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VdB  (0VdB=1  e-8  m/s)  VdB  (OVdBsl  e-8  nVs) 


1/1  Octave  Band  Analysis 
FF-1062  &  FF1071  Fire  Pump 
Each  Data  Set  Use  20  Averaged  FFT  Data 


(a)  Pickup  Location :  UMB 


(b)  Pickup  Location  :  LMB 


Figure  14  The  Summarized  1/1  Octave  Band  Levels  for  Type  III  Fire  Pumps 
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VdB  (OVdB=l  e-8  nVs)  VdB  (0VdB»1  e-8  m/s)  VdB  (OVdB-1*^  m/s) 


(c)  Pickup  Loctfion  :  LMB(A) 


Figure  14  The  Summarized  1/1  Octave  Band  Levels  for  Type  in  Fire  Pumps  (Cont) 
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Ikble  16:  SUMMARIZED  MEAN  (m)  AND  STANDARD  DEVUTION  (o)  OF  1/1 

OCTAVE  BAND  LEVELS 


Type 

Pickup 

Location 

Band  Level  (VdB) 

OVdBsle-8  m/sec 

1st 

2nd 

3rd 

4th 

Sth 

6th 

7th 

8th 

9th 

1 

MB(FE) 

D 

845 

75.5 

95.0 

IQI 

834 

81.3 

69.4 

703 

73.4 

O 

4.87 

3.59 

1.56 

1.09 

172 

3.34 

4.84 

4  16 

4.70 

MB(CE) 

Q 

789 

81.8 

95.5 

91.3 

IQJ 

72.5 

71.2 

75.1 

448 

5.49 

179 

154 

1.34 

142 

1.09 

1.95 

3.45 

PB(CE) 

D 

83.2 

859 

101  3 

984 

100.5 

990 

86.8 

834 

• 

3.42 

5.05 

2.90 

4.12 

168 

3.29 

154 

4.11 

4.56 

PB(FE) 

D 

850 

89  7 

1029 

95.3 

103.9 

963 

840 

827 

853 

a 

4  76 

6  14 

2.90 

317 

360 

3.96 

3.71 

489 

3.38 

PB(FE/A) 

a 

82  8 

86.3 

1016 

991 

1066 

1013 

892 

88.9 

86.3 

o 

341 

3.21 

4.41 

5.30 

365 

429 

267 

386 

348 

11 

MB(FE) 

a 

74.0 

75.0 

900 

89.5 

83.7 

83  2 

698 

688 

71.7 

a 

7.67 

1  99 

12  44 

1  37 

139 

3.73 

5.28 

4.93 

8.23 

MB(CE) 

a 

73.6 

77.3 

98.1 

928 

852 

862 

72.2 

708 

76.5 

o 

603 

048 

214 

4.01 

0.12 

0.60 

0.05 

0.57 

0.94 

MB(CE/A) 

5 

81.3 

81  9 

1029 

102.0 

898 

94.2 

820 

766 

5 

I  86 

3.44 

1106 

387 

1.71 

109 

3.09 

138 

1.72 

PB(CE) 

D 

864 

84.3 

1063 

97.6 

91.2 

98.4 

839 

763 

^23 

1.36 

3.72 

600 

5.37 

1.18 

119 

119 

159 

1.93 

PB(FE) 

D 

78.2 

81.2 

1086 

99  3 

913 

920 

85.9 

83.7 

82.1 

3.55 

4.13 

5.57 

4.88 

0.95 

1.26 

1.38 

135 

3.18 

III 

UMB 

D 

87.6 

79.1 

91.1 

89.2 

89.4 

84.3 

73.9 

77.4 

76.1 

3.67 

3.99 

5.29 

5.21 

3.70 

3.22 

496 

175 

6.76 

LMB 

n 

87.0 

81.7 

95.6 

94.0 

89.3 

85.7 

74.3 

79.7 

79.7 

4.19 

3.38 

3.42 

195 

3.72 

149 

2.79 

1.46 

3.18 

LMB(A) 

D 

87.1 

79.3 

94.7 

91.9 

87.2 

85.9 

74.4 

75.3 

80.4 

3.19 

3.35 

4.57 

171 

3.20 

2.93 

1.21 

4.06 

3.57 

UPB 

a 

91.4 

88.7 

106.3 

101  3 

99.3 

97.8 

85.2 

83.9 

4.55 

6.51 

6.57 

4.24 

8.18 

7.73 

7.86 

7.04 

7.98 

LPB 

D 

93.8 

93.3 

100.4 

100.6 

102.8 

99.1 

88.3 

83.2 

85.3 

-1. 

4.60 

466 

4.94 

121 

4.35 

4.83 

5.64 

6.99 

6.61 
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2.  Artificial  Fault  Siaulatioas 


a.  Fault  SiumtadoH 

In  OTder  to  assess  the  1/1  octave  band  method,  three  cases  of  artificial  fault 
simulation  were  performed.  Before  we  simulated  these  cases,  we  had  to  define  a  “fault”. 
For  simplicity,  we  assumed  a  fault  could  be  i^iproximated  with  S  bars,  as  shown  in  Figure 
IS.  The  bandwidth  of  the  fault  is  about  4.S  Hz. 


Figure  15  Fault  with  6VdB  Gain  at  Center  Bar 


Based  on  Eq  (8),  if  the  vibration  amplitude  doubles  with  a  linear  velocity 
scale  then  the  VdB  level  will  increase  6VdB  at  the  corresponding  frequency.  Therefore,  we 
defined  a  f:.ult  as  having  a  6VdB  gain  at  the  center  bar  and  use  four  bars  with  3  VdB  (half 
the  power  of  the  center  bar)  and  2VdB  at  the  side  bands  to  simulate  leakage  effects  around 
the  corresponding  frequency 
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b.  Simulation  of  A  Misalignment  Fault 

Misalignment  occurs  when  the  center  lines  of  two  shafts  are  <^set  or  meet  at 
an  angle.  The  characteristics  of  misalignment  in  a  spectrum  include: 

•  High  amplitude  axial  peaks  and  radial  peaks  at  1, 2  and  3  orders  of  shaft  RPM. 

•  Higher  harmonics  of  the  shaft  RPM  (greater  than  4  orders  of  shaft  RPM)  are 
generally  low  in  amplitude. 

Thus,  it  is  reasonable  to  add  the  6VdB  gain  fault  at  the  first  and  sectmd  (X^der. 
Figure  1 6(a)  and  16(b)  show  the  spectrums  before  and  after  imposing  the  fault.  Figure  16(c) 
compares  the  differences  in  the  1/1  octave  band  levels  and  the  10-5,000  Hz  broad  band 
level  for  these  two  conditions.  As  can  be  seen  in  Figure  16(c),  the  1/1  octave  band  levels 
have  a  4.46  VdB  gain  at  third  band  (which  corresponds  to  1  order  of  shaft  RPM)  and  a  3.13 
VdB  gain  at  the  forth  band  (which  corresponds  to  2  orders  of  shaft  RPM).  However,  the  10- 
5,000  Hz  broadband  level  only  increased  2.17  VdB.  Obviously,  the  1/1  octave  band  method 
is  more  sensitive  than  broadband  method. 
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VdB  (0VdB=1o-8m/s)  VdB 


Case  ^Misalignment  Fault 

CG-59  Fir©  Pump 
Pickup  Location :  PB(CE) 


(a)  Good  Condition  (b)  Damaged  Condition 


Orders  Orders 


(c)  1/1  Octave  Band  Levels 


io’  10*  10*  10* 

Frequency  (Hz) 


Figure  16  Simulation  Result  for  Artificial  Misalignment  At  Shaft 
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e.  Simulation  of  A  Looseness  Fault  At  Impdler 

Impeller  looseness  is  a  rotating  dement  looseness.  The  important 
characteristics  of  looseness  in  a  spectrum  include: 

1 .  Presence  of  a  large  number  of  hannonics. 

2.  Presence  of  half-barmonics. 

For  example,  the  fire  pump  of  a  CG-S9  has  six  impdlers  ctmtributing  to  a 
forcing  frequency  of  6  orders  of  the  shaft  RPM,  and  we  expect  higher  levels  at  6  orders  and 
its  harmonics  (12  orders,  24  orders,...  etc.)  fcx^  a  loose  impeller.  The  levds  at  half- 
harmonics  (e  g  3  orders  and  9  (vders)  will  also  increase.  Thus,  in  this  looseness  fault 
simulation,  an  artificial  fault  was  applied  by  adding  a  fault  with  6VdB  at  center  bar  at  die 
6th  order(i.e.  6  x  shaft  RPM)  and  adding  two  fault  with  3VdB  at  center  bar  at  the  3rd  <xder 
and  the  9th  order.  Figure  17(a)  and  17(b)  show  the  spectrum  of  good  and  damaged  fire 
pumps  Figure  17(c)  compares  the  differences  in  1/1  octave  band  levels  and  the  10-5,000 
Hz  broadband  level  for  these  two  conditions.  At  the  fifth  of  the  1/1  octave  band  levels,  there 
is  a  1.5  VdB  gain.  However,  the  10-5,000  Hz  broadband  level  only  increased  0.29  VdB. 
Again,  the  1/1  octave  band  method  is  more  sensitive  than  broadband  method. 
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(SAU8-e^=gpA0)  QPA 


Case  II:  Looseness  Fault  At  Impeller 

CG‘59  Fire  Pump 
Pickup  Location :  PB(FE) 


(a)  Good  Condition 


Orders 


0  5  10 

Orders 


(c)  1/1  Octave  Band  Levels 


Frequency  (Hz) 


Figure  17  Simulation  Result  for  Artificial  Looseness  At  Impeller 
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d.  Simulation  of  A  Bearing  Fault 

For  a  steady  state  condition,  some  periodic  signatures  exist  which  relate  to 
corresponding  bearing  faults,  called  bearing  frequencies.  These  bearing  frequencies  can  be 
found  in  the  VTAG.  In  this  section,  an  artificial  fault  imposed  at  the  motor  bearing  (free 
end)  of  a  CG-S9  fire  pump  has  been  simulated.  Hie  bearing  frequencies  of  this  bearing  are 
tabulated  in  Table  17. 


TABLE  17:  BEARING  FREC^UENCIES  OF  MRC310  BEARING 


Bearing  Frequency 


Train  order  of  rolling 
element 


Relative  rotation  order  of 
rotating  raceway 


Spin  order  of  rolling 
elements 


Irregularity  order  of 
rolling  element 


Irregularity  order  of 
rotating  raceway 


Irregularity  order  of 
stationary  raceway 


Symbol 


Order 


Suppose,  for  example,  that  there  is  a  wear  degradation  in  the  inner  raceway. 
Then  the  presence  of  the  fundamental  fj^  tone  with  harmonics  would  be  expected.  Figure 
18(a)  shows  the  spectrum  of  a  good  bearing.  The  6VdB  fault  is  imposed  at  the  4.9S  order 
as  shown  in  Figure  18(b).  Figure  18(c)  shows  the  resulting  increase  in  the  1/1  octave  band 
levels.  The  4.95  order  is  located  at  the  fifth  of  the  1/1  octave  bands.  The  10-5,000  Hz 
broadband  levels  before  and  after  damage  are  almost  unchanged.  However,  the  fifth  octave 
band  level  has  a  1 .57VdB  increase.  This  illustrates  why  a  single  broadband  level  cannot  be 
used  to  determine  the  condition  of  a  machine. 
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VdB  (0VdB=1a-8m/s) 


Case  lll:Bearing  Fault 

CG-59  Fire  Pump  (1 1/30/89) 
Pickup  Location :  MB(FE) 


(a)  Good  Condition 


Orders 


(b)  Damaged  Condition 


Orders 


(c)  1/1  Octave  Band  Levels 


Frequerwy  (Hz) 


Figure  18  Simulation  Result  for  Artificial  Bearing  Fault  At  Motor 
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V.  RESULTS  OF  USING  NAVSSES  FREQUENCY  DOMAIN 

DATABASE 


Since  digital  vibration  survey  instruments  are  more  popular  and  easy  to  use  than  other 
types,  only  frequency  spectrum  reading  were  recorded  and  stored  into  the  database.  Thus, 
our  1/1  octave  band  analysis  was  performed  using  frequency  domain  data.  Data  mea»ired 
for  the  CG-S9  fire  pumps  in  six  survey  dates  are  available  in  ASCII  format.  The  total 
number  of  fire  pumps  analyzed  was  twenty-seven.  Table  18  shows  the  available  frequency 
data  and  corresponding  survey  date.  Each  set  of  data  contains  a  10  order  and  a  100  order 
spectra  The  resolution  for  each  set  of  data  is  400  lines.  The  averaged  spectra  using  these 
data  are  shown  in  Figure  19  and  Figure  20  for  each  order  range.  The  thin  lined  curve 
represents  the  mean  levels  and  the  thick  lined  curve  represents  for  the  mean  plus  one 
standard  deviation  levels.  The  10  order  spectrum  was  used  to  compute  the  first  six  1/1 
octave  band  levels  and  the  100  orders  spectrum  was  used  to  compute  the  seventh  and  above 
band  levels.  Results  of  the  1/1  octave  band  analysis  fo  *  each  set  of  data  are  shown  in 
Appendix  F.  A  summary  of  1/1  octave  band  analysis  results  for  each  pickup  location  using 
all  the  fire  pumps  is  shown  in  Figure  21.  As  can  be  seen  from  Figure  21,  the  higher 
frequency  levels  (seventh  band  to  ninth  band)  are  about  95  VdB  for  the  motor  bearing  and 
1 1 0  VdB  for  the  pump  bearing.  However,  these  higher  frequency  levels  in  Figure  12  using 
time  waveform  tape  data  are  about  80  VdB  for  the  motor  bearing  and  90  VdB  for  the  pump 
bearing  As  we  reviewed  the  narrow  band  plot  for  the  NAVSSES  frequency  domain  data, 
we  found  that  the  levels  below  54  VdB  have  been  cut  off  and  forced  to  54  VdB.  This  is  why 
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levels  from  using  time  wtvefonn  ttpe  data.  Figure  22  diowt  a  food  cawnple  of  a  cut  off 
spectrum  for  a  motor  bearing  using  the  NAVSSES  frequency  domain  data. 


Ihble  18:  DATA  COLLECTION  SHEET  FOR  AVAILABLE  FREQUENCY 
DOMAIN  DATA  FROM  NAVSSES  DATABASE 


Survey  Date 

Number  of  Fire 
Pumps 

Number  of  Pickup 
Location 

Number  of  Dau 
Set 

12-MAY-93 

6 

s 

30 

23-FEB-92 

5 

5 

25 

21-OCT-90 

5 

5 

25 

07-JUN-90 

4 

5 

20 

26-FEB-90 

1 

5 

5 

30-OCT-89 

6 

5 

30 
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VdB  <0VdB=1 


Narrow  Band  Average  for  0-10  Order 
Using  NAVSSES  Database 
CG-59  Fire  Pump 
Each  Position  Use  27  Pumps 


(a)  Pickup  Location  ;  MB(FE) 


Figure  19  The  Averaged  Spectra  for  CG-59  Fire  Pumps  (0-10  Order)  Using 

Frequency  Domain  Data 
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VdB  (0VdB=1  e-8m/sec)  VdB  (0VdB=1  e-8m/sec)  VdB  (OVdBsl 


(c)  Pickup  Location :  PB(CE) 
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(d)  Pickup  Location  :  PB(FE) 


(«)  Pickup  Location :  PB(FE/A) 


Figure  19  The  Averaged  Spectra  for  CG-59  Fire  Pumps  (0-10  Order)  Using 

Frequency  Domain  Data  (Cont.) 
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(b)  Pickup  Location :  MB(CE) 
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Figure  20  The  Averaged  Spectra  for  CG-59  Fire  Pumps  (0>100  Order)  Using 

Frequency  Domain  Data 


(a)  Pickup  Location :  MB(FE) 
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VdB  (0VdB=1  e-8m/sec)  VdB  (OVdBsl  a-8m/sec)  VdB  (OVdB«1 


Figure  20  The  Averaged  Spectra  for  CG>59  Fire  Pumps  (0-100  Order)  Using 

Frequency  Domain  Data  (Cont.) 
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VdB  (0VdB=1e-8  nVs)  VdB  (0VdB=1  a-8  m/s) 


1/1  Octave  Band  Analysis 
Using  NAVSSES  Database 
CG-59  Fire  Pump 
Each  Position  Use  27  Pumps 


(a)  Pickup  Location  :  MB(FE) 


Frequency  (Hz) 


(b)  Pickup  Location  ;  MB(CE) 


Frequency  (Hz) 


Figure  21  The  Summarized  1/1  Octave  Band  Levels  for  CG>59  Fire  Pumps  Using 

NAVSSES  Frequency  Domain  Data 
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VdB(0VdB=1e-8nfVs)  VdB  (0VdB=1e-8  m/s)  VdB(0VdBs1e-8m/s) 
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(d)  Pickup  Location  :  PB(FE) 
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Figure  21  The  Summarized  1/1  Octave  Band  Levels  for  CG-S9  Fire  Pumps  Using 
NAVSSES  Frequency  Domain  Data  (Cont.) 


Frequency  (Hz) 


(e)  Pickup  Location  :  PB(FE/A) 


(c)  Pickup  Location  ;  PB(CE) 
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VdB  (0VdB=1«-8m/sec)  VdB  (0VdB=1a-8tn/ssc)  VdB  (0VdBs1e-8m/s«c)  VdB  (0VdBsl«-6m^sac) 


(a)  0-10  ordar  (or  MB(FE) 
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Figure  22  The  Narrow  Band  Spectra  for  CG-59  Fire  Pumps  Using  NAVSSES 
Frequency  Domain  Data  Measured  at  MB(FE)  and  MB(CE) 
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(d)  0-100  order  for  mB(CE) 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


Both  time  and  frequency  domain  analyses  were  performed  using  fire  pump  vibration 
data  to  determine  an  appropriate  alert  level.  Cr(»s  Peak  Analysis  in  the  time  domain  and  1/ 
1  Octave  Band  Analysis  in  the  frequency  domain  were  used  to  set  die  alert  threshold 
level(s).  The  following  conclusions  are  drawn; 

Time  Domain  Analysis  (Cross  Peak  Analysb) 

a.  The  measured  peak  envelop  data  at  five  different  pickup  locations  follows  a 
Gaussian  probability  distribution  in  the  VdBOog)  domain  well. 

b.  The  |i+2a  value  computed  using  a  peak  envelop  probability  density  function  in  the 
VdB  domain  gives  a  broadband  peak  amplitude  alert  level  rather  than  the  energy  content  in 
the  vibration  signal. 

c.  We  used  12  measured  data  sets  (6  fire  pumps  and  2  different  measurement  dates) 
which  were  available  on  cassette  data  tapes.  To  improve  the  quality  of  averaging,  more 
measured  data  sets  would  be  required. 

d.  NAVSSES  does  not  save  time  domain  data  on  any  storage  media  at  this  time. 
However,  for  some  shipboard  machinery,  such  as  low  and  high  pressure  air  compressors, 
time  domain  data  is  required  to  cast  measured  data  in  the  time-frequency  domain  to  detect 
the  possible  faults. 

Frequency  Domain  Analysis  (1/1  Octave  Band  Analysis) 

a.  1/1  Octave  Band  Analysis  (OBA)  uses  nine  frequency  band  alert  levels  which 
provides  more  detailed  information  about  machine  condition  than  the  simple  broadband 
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frequency  band  or  bands,  narrow  band  zoom  mode  analyas  can  be  perfcMmed  for  the 
selected  frequency  band(s)  to  identify  the  c(mip^ent(s)  which  may  have  faults 

b.  OB  A  divides  the  frequency  range  into  10  bins  over  lOkHz.  The  VdB  level  in  each 
frequency  bin  is  quite  sensitive  to  changes  in  the  energy  content  of  the  measured  vibratim 
signal. 

c.  The  OBA  approach  can  be  easily  incorporated  into  the  Digital  Vibratitm  Survey 
Instrument  (DVSI). 

d.  Five  sets  (27  fire  pumps)  of  400  line  frequency  domain  DVSI  data  were  received 
from  NAVSSES  and  OBA  was  performed.  The  VdB  levels  in  the  higher  frequency  bins 
were  generally  high  This  was  caused  by  the  fact  that  the  vibration  signals  which  were 
lower  than  the  VdB  limit(S4  VdB)  were  set  to  the  lower  VdB  limit  in  a  specified  dynamic 
range. 

e.  OBA  was  also  performed  using  time  domain  data  obtained  from  cassette  data  tapes. 
The  changes  in  VdB  levels  at  a  particular  frequency  and  the  corresponding  sidebands  were 
introduced  in  narrow  band  FFT.  The  results  are  transferred  into  the  octave  band  domain  and 
the  sensitivity  in  VdB  levels  in  frequency  bins  was  analyzed.  The  results  are  quite 
promising. 

f  NAVSSES  maintains  400  line  frequency  domain  data  for  all  machineries  in  a 
graphics  postscript  format,  not  in  ASCII  or  digitized  format.  It  was  not  possible  to  convert 
this  graphics  format  data  to  ASCII  format. 

Based  on  current  studies,  the  following  recommendations  are  made; 
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saved  using  DVSI.  However,  for  some  shipboard  madiinefy  sudi  as  low  and  lu^  pressure 
air  compressm-s,  time  domain  data  is  required  to  cast  the  measured  data  into  the  time- 
frequency  dmnain  to  detect  the  possible  faults.  We  suggest  saving  time  dcmiain  data  fw  the 
limited  number  of  shipboard  machinery. 

b.  We  suggest  implementing  the  1/1  Octave  Band  Analysis  Method  in  DVSI. 

c.  NAVSSES  maintains  400  line  frequency  domain  data  for  all  machinery  in  graphics 
postscript  format,  not  in  ASCII  or  digitized  format.  It  was  not  possible  to  cmivert  this 
graphics  format  data  to  ASCII  format  We  suggest  saving  future  files  in  ASCII  format  to 
permit  further  analysis  of  the  data. 

d.  We  suggest  saving/storing  DVSI  data  “as  is”,  without  alteration. 
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APPENDIX  A.  MATLAB  PROGRAM  CODE 


A.  ON-LINE  HELP  DOCUMENTATION 

This  section  of  the  guide  contains  all  of  the  MATLAB  functions  that  we  have 
developed  during  this  research.  Each  function  contains  a  purpose,  synopsis  and  description 
in  alphabetical  order.  The  MATLAB  built-in  functions  used  in  these  programs  are 
described  in  the  MATLAB  reference  guide. 

The  code  was  written  for  the  MATl  AB  4.0  version.  For  those  who  want  to  use  these 
programs  in  older  versions  need  to  consult  the  MATLAB  reference  guide  for  necessary 
modifications. 
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FiiriMM 


kurtosis  or  fourth  momoit  vectors  tnd  matrices. 

SyBQPsia 

kurtosis=kurt(x) 

Description 

krut  calculates  the  kurtosis  or  fouth  moment  value. 

For  vectors,  kurt(x)  is  the  kurtosis  value  of  the  dements  in  vecttx’  x.  For 
matrices,  kurt(x)  is  a  row  vector  containing  the  kurtosis  value  of  each 
column. 
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Purpose 


Plot  I /I  octave  band  bar  plot 

Synopsis 

logbar(Vdb) 

Iogbar(Vdb,LineStyIe) 

logbar(Vdb,LineStyle,UneWidth) 

Description 

logbar(Vdb)  draws  a  1/1  octave  band  bar  plot  using  Vdb  as  the  y  axis 
and  the  ANSI  prefetred  center  frequencies  and  pass  bands  as  the  x  axis. 
LineStyle  and  LineWidth  are  optinal  parameters.  The  default  LineStyle 
is  a  solid  line  See  “plot”  for  more  information  about  LineStyle  The 
default  LineWidth  is  O.S. 


% 


65 


Statistical  analysis  for  time  domain  data. 

Synopau 

mstat(y,fs,nbar,id) 

mstat(y.fs,nbar,id.manaly,mid,locat,mregen,pri) 

PtscriptiQB 

mstat(y,fs,nbar,mid,locat,manaly,mthre,z,mregen,pn)  will  analyze  the 
original  data^  using  the  manaly  method  and  produce  a  statistical 
report  with  a  nbar  pdf  plot. 

mstat(y)  allows  the  user  to  execute  the  program  in  an  interactive  way. 
y  amplitude  (calibrated  velocity  data) 

fs  sampling  frequency  (Hz) 

nbar  number  of  bars  in  Probability.  Density  Distribution  Plot 

mid  type  of  machinery 

locat  pickup  location 

manaly  method  of  analysis 

means  All  data  points  analysis 
-2  means  Overall  peak  points  analysis 
=3  means  Cross  peak  points  analysis 
mthre  with  threshold  or  not 

~\  without  threshold 
=2  with  threshold 
z  percentage  of  threshold 

mregen  transform  to  Vdb  or  not 

=1  without  transform 
=2  transform  to  Vdb  domain 
pri  print  or  not  (‘y’es  or  ‘n’o) 
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Furpoac 


Cross  peak  points. 

Synopsis 

ypc=npc(y) 

Description 

ycp=npc(y)  will  find  the  peak  points  between  each  two  zero-crossing., 
y  original  data 

yep  cross  peak  data 
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Overall  peak  points 

Synopafa 

yp^peak(y) 

PMcription 

yp=mpeak(y)  will  find  the  ov^-all  peak  points, 
y  original  data 

yp  overall  peak  data 


fttrofltf 


Threshold 

Sjjumsis 

yt=nthre(y,z) 

Description 


yt»nthre(y,z)  will  remove  the  points  with  z%  threshdd  . 
y  original  data 

yt  data  after  threshold 

z  percentage  of  the  direshold 

If  z=l  then  the  program  will  remove  the  peak  data  small  than 
1%  of  the  maxinum  peak  value. 
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1/1  octave  band  level. 


SXBIIPSil 

oct*octave(y) 

oct*octave(y, 11,12) 

Dcacriptiaa 

oct*(XTAVE(y)  returns  the  1/1  octave  band  Vdb  values  from  band  1 
to  band  10. 

oct=OCTAVE(y, 11,12)  returns  the  1/1  octave  band  Vdb  values  from 
band  i  i  to  band  12  and  the  others  bands  will  be  assigned  a  value  of  0. 
y(:,l) ;  Vdb  values 
y(:.2) ;  frequencies 


Purpose 


To  evaluate  the  avoaging  VdB  values  using  banning  window. 

Sjjumsis 

VdB=‘psd(y,nblock) 

PttcriptiQD 

VdB»psd(v,n&/ocit)  returns  the  average  VdB  values  using  time  domain 
data  y.  The  number  of  avera^ngs  is  q)eciried  by  nblock.  The 

recommend  number  of  averagings  is  20.  The  data  matrix  y  must  be  a  2° 
by  2  matrix.  41  contains  time  series  tnd y(:.2J  contains  conreqxmding 
velocity  dau.  The  banning  window  is  used  to  damp  out  the  dfects  of  the 

Gibbs  phenomenon.  The  r^erence  OVdB  is  10*^  m/sec. 
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Skewness  or  third  moment 


Sypopau 

skewness’=dcew(x) 

PMcrintion 

skew  calculates  the  skewness  ox  third  moment  value. 

For  vectors,  skew(x)  is  the  skewness  value  of  the  elements  in  vecux  x. 
For  matrices,  skew(x)  is  a  row  vector  containing  the  skewness  value  of 
each  column 


Purpose 


Plot  time  series  and  frequency  spectrum. 

SyoQpsis 

[Ya,y,fx]=transfer(yo,t,fo,nblock,sec) 

Description 

transfer(yo,t,fo,nblock,sec)  return  three  vectors,  plot  time  series  and 
frequency  spectrum.  The  return  vectors  are: 

Ya  averaged  fft  vectm*  of  yo 

y  calibrated  velocity  vector 

fx  frequency  vector 

The  input  parameters  include: 
yo  original  amplitude 

t  time  series  corresponds  to  yo 

fo  operating  frequency  (Hz) 

nblock  number  of  data  blocks 
sec  interested  frequency  sections 

The  spectrum  will  indicate  the  peak  position  assign  in  sec. 
For  example,  input  sec=(5  150],  the  program  will  search  the  peak  point 
between  5  Hz  and  150  Hz. 
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1.  kiiit.ai 

ftinctioo  kurtosis  -  kuft(y) 

SKURTKiulosis. 

%  For  vectors.  KURT(y)  ictums  the  kurtosis. 

%  For  mstriocs,  KURT(y)  is  a  row  vector  containing 
%  the  kurtosis  of  each  column.  See  also  COV 

%  Chao-Shih  Liu  4-26-93 

(n.m]=si2e(y); 
if  n»*l 

y=y’; 

{n.m]»size(y); 

end 

for  i“l:m 

meany  (;.i)=mean(y(;.i)).  ♦ones(n.  1 ). 
stdy  (:,i)»^d(y(:.i))*ones(n,  1 ); 
end 

s*((y -meany ). /stdy  ).  M ; 
kurtosis  =  sum(s)/n, 
end 
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function  k>gbar(VdbJ.ineStyleJ.ineWidth) 

% 

%  logbar  Plot  I/l  octave  band  bar  plot  using  aemilog  in 
%  X  axis(freq.)  and  linear  scale  in  Y  axis(Vdb) 

% 

%  logbar(Vdb)  plots  the  semilog  bar  plot 
%  for  1/1  octave  band  analysis. 

% 

%  logbar(VdbXineStyle)  plots  the  semilog  bar  plot 
%  for  1/1  octave  band  analysis  with  various 
%  types.  See  “plot**  f<v  more  about  LineStyk 
% 

%  Line  Width  is  a  optional  parameter  to  control  the  line 
%  width  of  the  bar.  Default  value  is  0.5 

y. 


%  Written  By  Liu,Chao-Shih  7/27/93 

if  nargin<2 
lineStyle=*-’. 
end 

if  nargin=  1 
linewidth=0.5; 
end 


% 

%  Define  the  upper  and  lower  limit  of  the  pass  bands  (Hz) 
% 

now=[11.2  22.4  45  90  180  355710  1400  2800  5600], 
lhigh=[now(2:10)  10000], 


% 

%  Determine  the  coordinates  of  the  bars 
% 

xl=flow(;); 

xh*fhigh(;); 

y=Vdb(:); 

(n,m]=si2e(y); 

n3=n*3, 

yy=zeros(n3+l,m); 

xx=yy, 

xx(l;3:n3,:)=xl(l:n), 

xx(2:3;n3,-.)=xh(l  n), 

xx(3;3;n3,:)=xh(l:n), 

xx=(xl(l);xx], 

yy(l:3:n3,:)=y, 

yy  (2.3.  n3,  :)=>•; 

yy=(0;yy]; 

set(semilogx(xx,yy,LineStyle),’linewidth’J.ine  Width) 
end 
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ftncticD 

% 

%  insUt(y/s.iibtfj&idJocat.iiiaMly;mdue.Z4nKten4sh) 

H  will  aiuly  sis  Uw  origiiul  daU  y  and  produce  a 
%  statistical  report  wi^  a  nbtf  pd[  plot. 

% 

%  instat(y)  will  allow  the  user  to  peifonn  the  program  in  an 
%  interactive  way. 

% 

%  y  ;  amplitude  (after  calibration) 

%  fs  :  sampling  frequency  (Hz) 

%  nbar :  number  bars  in  Prob.  Density  Distribution  Plot 

%  mid  :  type  of  machinery 

%  locat ;  pickiq)  location 

%  manaly;  method  of  analysis 

%  *'1  means  All  data  points  analysis 

%  *2  means  Overall  peak  points  analysis 

%  =3  means  Cross  peak  points  analysis 

%  mthre  :  with  threshold  or  not 

V«  ~1  without  threshold 

Vo  =2  with  threshold 

%  z  :  percentage  of  threshold 

%  mregen:  transform  to  Vdb  or  not 

%  =1  without  transform 

%  =2  transform  to  Vdb  domain 

%  pri :  print  or  not 

Vo  ='y’  means  to  print  out  the  result 

% 

Vo  Liu,Chao-Shih  6/7/93 
Vo 


if  nargin==l 

fs»input(*  Input  the  sampling  frequency  ?  (Hz)’); 
nbar=input(‘Input  the  number  of  bar  used  in  the  PDF  plot  ?’), 
mid*input(‘ Input  the  type  of  the  machinery  ?’,’s’); 
locat-input(‘ Input  the  pickup  location  ?’,’s’), 
manaIy*menu(‘Choose  the  method  of  analysis',... 

‘All  data  points  analysis’,.  . 

‘Overall  peak  analysis',... 

‘Cross  peak  analysis’); 
mthre=menu(‘Type  of  the  threshold  ?’,... 

‘Without  threshold’,... 

‘With  nVo  threshold’); 
if  mthie*=2 

z»input(‘ Input  the  %  of  the  threshold’);z=z/100, 
end 

mregen=menu(‘Transfotm  the  data  ?’,... 

‘Original  data’,... 

‘Vdb  Domain’); 
end 

% 
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TV  1VSVU1IM  VI  simjaiv 
% 

if  nianaly=l 

int=’All  data  points’ 
elseif  inanaly==2 
yp=npeak(y), 
int“’Overall  peak’, 
else 

yp=npc(y); 
int= ’Cross  peak’; 
end 

% 

%  Type  of  threshold 
% 

if  mthre==l 

tht=’  without  threshold’; 
else 

tht=(‘  with  ‘  num2slr(z*100)  ’%  threshold’], 
yp=nthre(yp,z); 
end 

pnuniber=length(yp); 


% 

%  Compute  the  first  four  mements  for  linear  domain 
% 

yp=abs(yp); 

mepx=(‘m=’  num2str(mean(yp))); 
stdpx=[‘s=’  num2str(std(yp))], 
ske wpx= ( ‘  Skewness=  ’  num2str(skew(y p))] , 
kurtpx=  [  ‘Kurtosis=  ’  num2str(k>.irt(}’p))]; 

% 

%  Compute  the  first  four  moments  for  Vdb  Domain 
% 

if  mregen==2 
yr=20*logl0(yp/le-8); 
meany  r=mean(y  r);stdy  r=std(>'r); 
skew>'r=skew(yr)>curtyr=kurt(yr); 
merpx=['m  =  ‘  num2str(mean0’r))l, 
stdrpx=[‘s  =  ‘  num2str(std(yr))], 
skewrpx=( '  Skewness=  ’  num2str(skew(y r))], 
kurtrpx= [  ‘  Kurtosis=  ’  num2sU'(kurt(y  r))] ; 
end 


% 

%  Plot  the  overall  envelope  PDF 
% 
clg, 

if  mregen==l 
subplot(312) 

[av,au]=hist(yp,nbar);bar(au,av); 
y label(‘Probability  ’);xlabel(‘  Velocity  (m/sec)’); 
title(’Original  Data(x)  Envelope’), 
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■el(tex^0.3.1  .nepxX  YailB\’iygdb(d\t<)iM»’/i^*); 
tet(text(0.3.0.8jidpxX’foola\*iiymM\Yo^’.’italic‘); 
aet(text(0.3.0.6^ewpx).’foola\’tinies’.YooU’.’iUlic’); 
•et(texK0.3.0.4Jniitpx).Yoam\’tii]ies’.'fooU’.’iUlic'); 
else 

subplot(323) 

[av.au]«faist(yp^bv);bar(au,av); 
ylabelCPiobability  ’);xlabel(' Velocity  (m/sec)’); 
tiUe(‘(a)  Linear  D^ain*); 
subplot(324) 

[av.au]*hist(yr4)bar);bar(au,av); 
ylabelCPiobability  •);)dabel(‘ Vdb  (0Vdb=le-8  m/scc)’); 
title('(b)  Vdb  Domain’); 
subpIot(3 1 3);axis(‘off ’); 

set(text(0. 1 . 1  ^epxX’fontn’.’symbol’.’fonta’.’iUlic’); 
set(text(0.1,0.8,stdpx),’fontn’.’symbor.‘fonu*.‘italic’); 
sel(text(0. 1 .0.6,skewpx).  ’fontn  ’times’.’fonU  ‘.’iulic  ’). 
set(text(0. 1 .0.4  Jcurtpx),’fontn’.’times’,’fonta*,’italic’); 
sel(text(0.7,l.nierpx),’fontn’,’symbor,’fonta*,’iulic’), 
set(text(0.7,0.8,stdrpx). ’fontn’, ’symboi’.’fonu’, ’italic’); 
sel(text(0.7,0.6,skewTpx),’fontn’.’times’,’fonta’,’italic’); 
set(text(0.7,0.4Jcurtrpx),’fontn’,’Umes’.’fonta’.’italic’); 
end 


% 

%  Statistical  analysis  results 
Vo 

subplot(3 1 1  ),axis(‘off ’), 

tal=['StatisticaI  Analysis  Results  for  '  mid], 

set(text(0.5,l,tal),’fontsi’,14.’fontw’,’bold’,’hori’,’c’); 

loc=(‘ Pickup  Location  :  *  locat], 

set(text(0.5,0.8,loc).’fontsi’,14,’fontw’,‘bold’,’hori’,’c’); 

ta2=['Data  analysis  by  using  :  ’  mt  tht], 

set(texl(0.5,0.4,ta2),’hori’,’c’); 

block=[‘Number  of  data  points  :  ‘  num2str(pnumber)]; 

sel(texl(0.5,0.2,block),’hori’,’c’); 

frange=[‘The  sampling  frequency  :  ‘  num2str(fs)  ‘  Hz’]; 

set(tcxl(0. 5,0,frange),  ’hori  ’c  ’ ); 


% 

%  Print  the  graph  or  not 
% 

if  nargin<4 

pri=input(' Print  the  plot  ?  (y/n)’,’s’), 
end 

if  pri=’y’ 

set(l,’PaperPosition’,(1.25  2  6  8.5]) 
print  -fl 
end 
end 
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'WI“TipCIHr 


functicm  yq>^pc(y) 

% 

%  ycp^pc(y)  will  find  the  peak  points  between  each  two 
%  zero-crossing 
% 

%  y  .  original  data 
%  yep :  cross  peak  data 
% 

% 

%  Liu.Chao-Shih  4/26^3 
% 

y-y(:). 

k*°length(y); 

y  *y  -mean(y  )*ones(k,  1 ), 

npc“0;l*l; 

ycp=zeros^.l); 

for  n=2;k 

if  sign(y(n)*y(n-l))==-l 
h=n-l  ;npc=npc+l ; 

(y  p.np]=max(abs(y  0 :  h))); 
np=np+l-l, 
ycp(npc)=y(np); 
l=h+l; 
end 
end 

ycp=ycp(  1  :npc); 
end 
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f  mctioo  yp*iiqpeak(y) 

% 

H  yp^peak(y)  will  find  the  overall  peak  points. 
% 

%  y  ;  oiigina]  data 
%  yp :  overall  peak  data 
% 

% 

%  Liu.Chao-Shih  4/14/93 
% 

y-y(). 

k-length(y); 

y=y-mean(y)*ones(length(k),  1 ); 
yp^ros(k,l); 

yp(i)=y(i); 

pc=l; 

forn=2:k-l 

nsl=n-l; 

nl  =sign((y(n+l  )-y(n))*(y(n)-y(nsl ))); 

ifnl**-l  |y(n)*%(nsl) 
pc=pc+l;yp(pc)=y(n). 
end 
end 

yp=yp(l:pc,l). 

end 
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w. — Biuie.ni 


function  yt“nthre(y,z) 

% 

%  yt=nthre(y,z)  will  remove  the  points  with  threshold 
% 

%  y  :  original  data 

%  yt :  data  after  threshold 

%  z  :  percentage  of  the  threshold 

%  If  z~l  then  the  program  will  not  return  the  peak 

%  peak  data  small  thw  1%  of  the  maxinum  peak  value 

% 

%  Liu.Chao-Shih  4/14/93 
% 

y*sort(y(:)). 
ly*length^); 
maxy=max^); 
miny =z*0.0 1  •maxy . 
for  1=1  :ly 
if  yO)>=miny 
po=l; 

yt=y(po:ly.l); 

break 

end 

end 
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fiweHoo  oci«oGicve(y4l42) 

K 

%  OCTAVE  Octave  band. 

% 

%  octKXTAVE(y)  ictunis  the  1/1  octave  band  Vdb  values  from 
%  band  1  to  band  10. 

% 

%  oct*OCTAVE(y,i  142)  returns  the  1/1  octave  band  Vdb  values 
%  from  band  il  to  band  i2  and  the  others 
%  bands  will  be  assigned  to  0. 

% 

%  y(..l) :  Vdb  values 
%  y(:^) ;  frequencies 

% 

%  Chao-Shih  Liu  7/13/93 
% 

fc=[16  31.5  63  125  250  500  1000  2000  4000  8000], 
now=[l  1.2  22.4  45  90  180  355  710  1400  2800  5600], 
fhigh=[22.4  45  90  180  355  710  1400  2800  560011200], 
np=max(size(y)); 
oct=zcros(10,l). 


% 

Vo  Default  the  il  and  i2 
Vo 

if  nargin=  1 
il=l,i2=10; 
end 

Vo 

Vo  if  Vdb  lower  than  0  db,  force  the  value^Odb 
Vo 

for  k*  1  :np 
if  y(k,l)<0 
y(k,l)=0; 
end 
end 


% 

Vo  convert  to  linear  velocity  value 
Vo 

vref=le-8, 

y(..l)-l0.^(y(.,l)/20)*vref, 
finax=y  (np,2),df=fmax/(np'  1 ); 

% 

Vo  calculate  the  1/1  octave  band  levels 
Vo 

df2=df/2, 

fupper=y(;,2)+df2;  %  Upper  limit  of  frequency 
flower=*y(;,2)-df2;  %  Lower  limit  of  frequency 
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fork-il:i2; 

marie  I  *fi4^>er-flow(k). 
inarK2*flower>flow^); 
mark3«fupper-fhig})(k), 
mark4-fk>wer-fhigh^X 

forj*l:np 

if  maric40)>^  *•  larger  than  the  band(k) 

break 

else 

if  maiklOK^  %  smaller  than  the  band(k) 

%  do  nothing 

else  %  cover  the  band(k) 

if  mark3(j)>0 

if  maik2(j><0  %  band(k)  smaller  than  df 
suni(k)-y  0. 1  r2*(fhi^(k)-flow(k)). 
break; 

else  %  partial  band  in  band(k-»-l) 

8um(k)=sum(k)+y(j.  1  )''2*abs(mark4(j)); 

end 

else 

if  maik2(j)<0  %  paitial  band  in  band(k-l) 

sum(k)=y  0, 1  )"'2*maric  1  (j); 

else 

sum(k)=sum(k)+y(j.  1  )''2*df, 

end 

end 

end 

end 

end 

oct(k)=20*log  1 0(sqrt(sum(k))/vreO. 
end 

bar(oct) 

xlabc!(' Number  of  Center  Frequency’) 
ylabel(' Vdb  value  (0Vdb=le*8  m/sec)’) 
axis((il-l  i2+l  40  120]) 
fork=il:i2 

octext=sprintf(‘%5.  If ’,oct(k)), 
set(tcxt(k,oct(k)+2,octext),’fontw’,’bold’,’horiz’,’ccnter’) 
end 
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t.  pMLn 

functioB 

%  Vdb-PSD(y.iibiock) 

% 

%  This  fuactioa  will  compute  the  $ivenge  Vdb  vtiue  by 
%  using  banning  to  nblodt  data  set 
% 

%  y  ;  original  daU  matrix  (time.velocity) 

H  nblock  :  number  of  data  blocks 
H  Ya  ;  average  lit  Vdb  vector  of  y(0Vdb*le-8  m/secj 
%  fx  :  frequency  vector 
% 

%  Liu.Chao-Shih  12: 1 1PM  7/14/93 

% 

%  standardize  the  data  matrix 

n-size(y,2); 

ifn~*2 

y=y‘; 

end 


t*y(;.i);y=y(  .2). 

dt=t(2)-t(l). 

fs=l/dt. 

vTef=le-8, 

% 

%  Compute  the  fft  data 
% 

n=length(y  ynblock,n2=n/2 , 
nl=l; 

w^haimingCn); 

fornb=l:nblock 

nh*nb*n, 

Y(  1  ;n^b)=ffl(w.  •y  (nl  :nh,  1  ))*dt, 
nl*nh+l, 
end 

if  nblock~^l 

Y=(sum(abs(Y)’)/nblock)’; 

else 

Y=abs(Y); 

end 

Ya=20*logl0(2*Y/vref). 

f=fs/2*(0:n2)/n2; 

Ya(n2+2:n)=[], 

fx-f(2:n2+l)’; 

df=fx(2)-fx(l); 

fx=lO;fxl. 

Vdb=[Ya4‘x]; 
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fuDctijn  skev^Dcss  « skevi{y) 

SSKEW  Skewness.  For  vectors,  SKEW(y)  returns  the  skewness 
%  For  matrices,  STD(y)  is  a  row  vector  coolainmg 
%  the  skewness  of  each  column.  See  also  COV. 

%Chao-Shih  Liu  4-26-93 

[n,m]«size(y), 
if  n**l 

y“y*. 

[n,m]-size(y); 
for  i»l:m 

meany(:,i)>nnean(y(:j)).*ones(n,l); 
sUty(:,i)»std(y(;4))*ones(n.  1 ), 
end 

s*((y  -meany ).  /stdy ).  ''B , 
skewness  =  sum(syn; 
end 
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ftuieiicii 

H  pf«j^]"«iMufciQro,t/o4iblodc;He) 

% 

H  This  functioo  will  plot  time  A  ficquency  dnmsin 
% 

%  yo  ;  origiiul  velocity  vector  (befoK  celibratiao) 

%  t :  time  vector 

%  fs ;  sample  frequenqr 

%  fo  ;  operation  frequency  ■  RPM/60 

%  nblodc :  number  of  data  blocks 

%  Ya  :  average  fit  ve^or  of  y 

%  lx  ;  frequency  vector 

%  y  :  original  velocity  vector  (after  calibration) 

V* 

%  Liu,Chao-Shih  I  ;S8  7/10/93 
% 

dt=t(2)-t(l); 

fs=l/dl; 

itext*input(‘ Input  the  title  of  this  record  =>’.V); 
iocat=input(input  the  location  of  dau  coUect’.’s’); 
y=yo; 

scale=menu( ‘Choose  the  method  of  X  scale’.’Frequency  (Hz)’,’Order’); 
vref=le-8; 


% 

%  Standarized  the  size  of  y 
% 

(xd,yd]=size(y), 
if  yd~=l 
y=y’;t=t'; 
end 

ny=length(y); 

% 

%  Compute  the  ITt  data 
% 

n=ny /nblock;n2=n/2 , 
nl=l, 

win=menu(‘Typc  of  Window  System ’.’Without  window  system’,  . 
‘Harming  window ’.’Hamming  w  indow’); 
if  win=l 
w=ones(n,  1 ); 
elscif  win=2 
w=hanning(n); 
else 

w=hamming(n); 

end 

for  nb=l;nblock 
nh=nb*n. 

Y(  1  :n.nb)=fft(w.  •y  (nl  :nh,  1  ))*dt, 

nl=nh+l, 

end 
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else 

Y-eb»(Y). 

end 

Ye-20*logl0(2*Y/vref). 

M‘s/2*(0;n2>h2; 

Ye(n2+2:n)-a; 

if  scek*~l 

fx-f(2:n2+l); 

df*fx(2)-fx(l). 

xt«’Frequency  (Hz)’;unit«’  Hz 
else 

fi(-f(2;n2+l)/fo; 

df-fx(2)-fx(l). 

k>w«ceil(0.7S/df)^igb*ceil(  1 .25/df). 
(Yinsx,Ypos]=in«x(YaOow.liigh));fscsle*fx(low+Ypos- 1 ); 
fx-fx^sc*le;fo^o*fsc*le, 

xt*[‘Frequency  (Normalized  Orders)’l;unit**  Order 
end 

df-fx(2Kx(l); 
sec”ceii(scc/df). 
if  sec(l)==<) 
scc(l)=l; 
end 

for  mk=l  :length(sec)-l 

[mav(mk),pos(mk)]=max(  Ya(sec(mk):  sec(mk+l ))), 
pos(mk)=pos(nik)+sec(mk); 
end 


% 

%  Plot  time  domain  and  frequency  domain 
% 

figuref  1  );clf,subplot(2 1 2), 

tf=I024; 

if  ny<1024 

tf*ny; 

end 

plot(t(l  :tf.  1  ),y(I  :tf,  1  ),’r’);axis((t(l  ,I )  t(tf,l )  min(y)  max(y))); 
titleCTime  Domain ’);xlabel(' Time  (sec)’);ylabel(‘Velocity(m/sec)’);grid 
% 

%  General  information 
% 

subplot(2 1 1  );axis(‘  off’ ) 

mfo=['The  operation  condition  =  ‘  num2str(fo*60)... 

‘  RPM  (=  ‘  num2str(fo)  ‘  Hz)’]; 
mloc=[‘Data  pickup  location  =  *  locat]; 
mfs=[‘Data  sampling  frequency  =  ‘  num2str(fs)  ‘  Hz’], 
mdp=[’Total  number  of  sampling  data  points  =  ‘  num2strGength(y))]; 
set(text(0.5,l.itext),’fontsi’,15,’fontw’,’bold’,’horiz’.’center’) 
text(0.05,0.9,mloc);text(0.05,0.824nfo); 
text(0.05,0.74,mfs),lcxt(0.05,0.66,mdp); 

text(0,0.5,’*  The  time  domain’); 

mt=[‘Mean  value  =  ‘  num2str(mean(y))  ‘  m/sec’]; 
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te]a(0.0S.0.4jitfX 

•^[‘Smidard  d^alioii  -  '  iium2str(ttd(y))  ‘  mAnc’); 
text(0.0S,0.32^); 

text(0.0S,0.24,(‘TlK  plot  only  shown  ‘  Duin2sti(tf)  *  data  pnints']), 
pauae 

nguic(2);clf;subplot(2 1 2); 
plot(fx.Ya(2:n2+l).’r’); 
iippef*(oeU((max(mav))/1 0)+l  )•  10. 
axis([0  lx(n2)  20  upper]). 

title(‘Frequency  Dmain’);xlabel(xt);ylabel(‘Velocity  in  Vdb').grid 

for  lop»l  :leng^(aec)-l 

text(fx(posOop)- 1  ).mav0op).num2sir(lop)). 

end 

% 

%  General  information 
% 

subplot(21  l);axis(‘ofr’) 

mfo=(‘The  operation  condition  =  ‘  num2str(fo*60)... 

‘  RPM  (“  ‘  num2str(fo)  ‘  Hz)’]; 
mloc*['Data  pickup  location  =  ‘  locat]. 
mfs=('Data  sampling  frequency  »  ‘  num2str(fs)  ‘  Hz’]; 
mbk=(‘Number  of  data  block  avaages  “  ‘  num2str(nbiock)]; 
mdp={‘Block  size  =  ‘  num2str(length(y)/nblock)  ‘  points/block’]; 
set(texl(0.5, 1  ,itext),’fontsi ’,1 5,’fontw’,’bold’.’horiz’.’center’), 
text(0.05,0.9,mloc);text(0  05,0.82.mfo), 
text(0  05,0.74 jnfs);text(0.05,0.66.mbk);text(0.05,0.58,mdp). 
text(0,0.5,’*  The  frequency  domain’); 
for  lop=  1 :  length(sec)- 1 

as=[num2str(lop),’  =  (*,num2str(fx(pos(lop)-l)),... 
unit,’,  *.num2str(mav(lop)),’  db)’], 
if  lop>4 

xpos=0.55;ypos=0.8-0  08*lop, 
else 

xpos=0.05,ypos=0.48-0.08*lop, 

end 

text(xpos,ypos,as), 

end 

morder=[‘l  order  =  ‘  num2str(fo)  ‘  Hz’], 
text(0.05,0.084norder); 

mref=[‘The  reference  0  Vdb  =  ‘  num2str(vref)  ‘  m/sec’]; 
te\t(0.05,0jnref), 

% 

%  zooming  up  the  plot  ? 

% 

zoom=l; 
while  zoom~=3 
subplot(212) 

zoom=menu(‘zoom  up  ?’,’0-10  order’,... 

‘Use  mouse  to  select  the  range’,’No’); 

ifzoom*=l 

axis([0  10  20  upper]); 

elseiif  2oom==2 

axis([0  fx(n2)  20  upper]); 
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xrange>sort(xs), 
if  xrajige(l>^ 
xrange(l)^, 
end 

•xis((xrange(l)  xrange(2)  20  un’O']). 

else 

break 

end 

end 

% 

%  Print  out  the  plot 
% 

uiiMTnt»(  *set(l,”PaperPosition”,[0.25  1  8  9.25J);*  ‘print  -fl;’ ... 
‘aet(2.”PaperPosition”.[0.25  1  8  9.251);*  ‘pnnt  -f2;*  ‘end’); 
uicontroi(‘riyle’,’push’.’iuiits*,’nonnar.’pos‘.(0.9  0.9  0.08  0.06], ... 
‘string ’[mnt*. ’call’, uiprint), 

uicontrol(‘style’,’push', 'units’, ’nom>at’.’pos‘.[0.9  0.82  0.08  0.06], ... 
‘string’,’exit’.’cair,’end‘). 
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APPENDIX  B.  EXAMPLES  OF  USING  MATLAB  PROGRAM  IN 

PC486 


MATLAB  allows  the  user  to  write  a  script  file  A  script  file  is  an  external  file  that 
contains  a  sequence  of  MATLAB  statements  The  following  example  will  show  the  script 
files  used  to  perform  the  analysis 

A.  CROSS  PEAK  ANALYSIS  EXAMPLE 

In  this  example,  a  CG-S9  fire  pump  measured  at  MB(FE)  data  set  is  used  in  the 
analysis.  The  file  name  of  the  data  file  is  flpl.asc  The  MATLAB  command  is  anchored 
by  ». 

» load  flpl.asc 
»  y=flpl(:,2), 

» t=flpl(.,l), 

»  mstat(y) 

The  user  can  then  interact  with  the  program  to  select  the  method  of  analysis.  Instead 

of  this  interactive  process,  the  user  can  use  the  following  to  execute  the  program  in  an 

automatic  way. 

» load  flpl  asc 
»  y=flpl(;,2), 
t=flpl(;,l), 

»  ypc=npc(y); 

»  yt=nthre(ypc,l); 

»  mid=’CG-59  Fire  Pump’, 

» Iocat=’MB(FE)’; 
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B.  1/1  OCTAVE  BAND  ANALYSIS  EXAMPLE 

Here,  proceeding  as  in  the  above  example,  a  CG-S9  fire  pump  measured  at  MB(FE) 
data  set  is  used  to  perform  the  1/1  octave  band  analysis.  The  file  name  of  the  data  file  is 
flpi.asc.  The  MATLAB  ccHnmand  is  anchored  by  ». 

» load  flpi  .asc 

»  yh*flpl(:,2);  %  yh  is  a  time  series  sampled  at  10  kHz. 

»  yl=yh(l:5:length(yh)),  %yl  is  a  time  series  sampled  at  2  kHz. 

» t=flpl(:,l); 

»  VdBh=psd(yh,20), 

»  VdBl=psd(yl,20); 

»  octh=octave(VdBh,7,9); 

»  octl=octave(Vdbl,l,6), 

»  octband=[octh  octl]; 

» logbar(octband); 
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APPENDIX  C.  FIGURES  OF  PROBABILITY  DISTRIBUTION  OF 

CROSS  PEAK  POINTS 

This  section  contains  all  the  cross  peak  distribution  figures. 
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Probability  Probability 


Probability  Density  Function  for  Cross  Peak  Points 
CG-59  Firs  Pump  (10/3(V89) 

Pickup  Location :  MB(FE) 


(a)  Fire  Pump  #1 


12  3  4 

Velocity  (m/sec) 


(b)  Fire  Pump  #2 


12  3  4 

Velocity  (m/sec) 


(c)  Fire  Pump  #3 


12  3  4 

Velocity  (m/sec) 


(d)  Fire  Pump  #4 


12  3  4 

Velocity  (m/sec) 


(e)  Fire  Pump  #5 


12  3  4 

Velocity  (nVsec) 


(f)  Fire  Pump  #6 


12  3  4 

Velocity  (nrVsec) 


Figure  Cl  The  Probability  Distributions  of  Cross  Peak  Data  for  CG-59( 10/30/89) 

Fire  Pumps  Measured  at  MB(FE) 
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Probability  Probability  Probability 


Probability  Density  Function  for  Cross  Peak  Points 
CG-59  Fire  Pump  (2/26/90) 

Pickup  Location :  MB(FE) 


(a)  Fire  Pump  #1 


(b)  Fire  Pump  «2 


12  3  4 

Velocity  (m/sec) 


(c)  Fire  Pump  #3 


12  3  4 

Velocity  (m/sec) 


(d)  Fire  Pump  #4 


12  3  4 

Velocity  (m/sec) 


(e)  Fire  Pump  #5 


12  3  4 

Velocity  (nVsec) 


(f)  Fire  Pump  #6 


12  3  4 

Velocity  (m/sec) 


12  3  4 

Velocity  (m/sec) 


Figure  C.2  The  Probability  Distributions  of  Cross  Peak  Data  for  CG-59  (2/26/90)  Fire 

Pumps  Measured  at  MB(F£) 
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Probability  Probability  Probability 


Probability  DanaKy  Function  for  Croaa  Peak  Pointa 
CG-59  Fira  Pump  (1(V3()l/89) 

Pickup  Location :  MB(CE) 
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Figure  C.3  The  Probability  Distributions  of  Cross  Peak  Data  for  CG>59(10/30/89) 

Fire  Pumps  Measured  at  MB(CE) 
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Probability  Probability  ProbabMity 


Probability  Density  Function  for  Cross  Peak  Points 
CG-59  Fire  Pump  (2/26/90) 

Pickup  Location :  MB(CE) 
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Figure  C.4  The  Probability  Distributions  of  Cross  Peak  Data  for  CG>S9  (2/26/90)  Fire 

Pumps  Measured  at  MB(CE) 
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Probability  Probability  Probability 


Probability  Density  Function  for  Cross  Peak  Points 
CG*59  Fire  Pump  (10/30/89) 

Pickup  Location :  PB(CE) 
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Figure  C.5  The  Probability  Distributions  of  Cross  Peak  Data  for  CG-59(1 0/30/89) 

Fire  Pumps  Measured  at  PB(CE) 
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Probability  Probability 


Probability  Density  Function  for  Cross  Peak  Points 
CG-59  Fire  Pump  (2/26/90) 

Pickup  Location :  PB(CE) 
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Figure  C.6  The  Probability  Distributions  of  Cross  Peak  Data  for  CG-59  (2/26/90)  Fire 

Pumps  Measured  at  PB(CE) 
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Probability  Probability  Probability 


Probability  Density  Function  for  Cross  Peak  Points 
CG-59  Firs  Pump  (10/30/89) 

Pickup  Location :  PB(FE) 
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Figure  C.7  The  Probability  Distributions  of  Cross  Peak  Data  for  CG-S9( 10/30/89) 

Fire  Pumps  Measured  at  PB(FE) 
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Probability  Probability  ProbabiHly 


Probability  Density  Function  for  Cross  Peak  Points 
CG-59  Fire  Pump  (2/26/90) 

Pickup  Location :  PB(FE) 
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Figure  C.8  The  Probability  Distributions  of  Cross  Peak  Data  for  CG-S9  (2/26/90)  Fire 

Pumps  Measured  at  PB(FE) 
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Probability  Probability  Probability 


Probability  Density  Function  for  Cross  Peak  Points 
CG-59  Fire  Pump  (10/30/89) 

Pickup  Location :  PB(FE/A) 
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Figure  C.9  The  Probability  Distributions  of  Cross  Peak  Data  for  CG>59(10/30/89) 

Fire  Pumps  Measured  at  PB(F£/A) 
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Probability  Probability  Probability 


Probability  Density  Function  for  Cross  Peak  Points 
CG-59  Fire  Pump  (2/26/90) 

Pickup  Location :  PB(FE/A) 
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Figure  CIO  The  Probability  Distributions  of  Cross  Peak  Data  for  CG-59  (2/26/90) 

Fire  Pumps  Measured  at  PB(FE/A) 
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Probability  Probability 


Probability  Density  Function  for  Cross  Peak  Points 

DDG-994  Fire  Pump 
Pickup  Location :  MB(FE) 


(a)  Fire  Pump  #1  (b)  Fire  Pump  #2 
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Figure  C.ll  The  Probability  Distributions  of  Cross  Peak  Data  for  DDG*994  Fire 

Pumps  Measured  at  MB(FE) 
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Probabiltty  Probability 


Probability  Density  Function  for  Cross  Peak  Points 

DDG-994  Fire  Pump 
Pickup  Location :  MB(CE) 


(a)  Fire  Pump  #1  (b)  Fire  Pump  #2 
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Figure  C.12  The  Probability  Distributions  of  Cross  Peak  Data  for  DDG-994  Fire 

Pumps  Measured  at  MB(CE) 
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Probability  Probability 


Probability  Density  Function  for  Cross  Peak  Points 

DDG-994  Fire  Pump 
Pickup  Location :  MB(CE/A) 
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Figure  C.13  The  Probability  Distributions  of  Cross  Peak  Data  for  DDG-994 
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Probability  Probability 


Probability  Density  Function  for  Cross  Peak  Points 

DDG-994  Fire  Pump 
Pickup  Location :  PB(CE) 


(a)  Fire  Pump  #1  (b)  Fire  Pump  #2 


(c)  Fire  Pump  #3 


Figure  C.14  The  Probability  Distributions  of  Cross  Peak  Data  for  DDG>994  Fire 

Pumps  Measured  at  PB(CE) 
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Probability  Probability 


Probability  Density  Function  for  Cross  Peak  Points 

DDG-994  Fire  Pump 
Pickup  Location :  PB(FE) 


(a)  Fire  Pump  #1  (b)  Rre  Pump  #2 
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Figure  C.15  The  Probability  Distributions  of  Cross  Peak  Data  for  DDG>994  Fire 

Pumps  Measured  at  PB(FE) 
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Probability  Probability 


Probability  Density  Function  for  Cross  peak  data 
FF-1062  Fire  Pump  (11/24/89) 

Pickup  Location :  UMB 
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Figure  C.16  The  Probability  Distributions  of  Cross  Peak  Data  fo.  FF*1062  Fire 

Pumps  Measured  at  UMB 
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Probability  Probability 


Probability  Density  Function  for  Cross  peak  data 
FF-1071  Fire  Pump  (3/10/89) 

Pickup  Location :  UMB 


(a)  Fire  Pump  #1  (b)  Fire  Pump  #2 


(c)  Fire  Pump  #3  (d)  Fire  Pump  #4 


Figure  C.17  The  Probability  Distributions  of  Cross  Peak  Data  for  FF-1071  Fire 

Pumps  Measured  at  UMB 
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Probability  Probability 


Probability  Density  Function  for  Cross  peak  data 
FF-1062  Fire  Pump  (11/24/89) 

Pickup  Location :  LMB 


(a)  Fire  Pump  #1 
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Figure  C.18  The  Probability  Distributions  of  Cross  Peak  Data  for  FF-1062  Fire 

Pumps  Measured  at  LMB 
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Probability  Probability 


Probability  Density  Function  for  Cross  peak  data 
FF-1071  Fire  Pump  (3/10/89) 

Pickup  Location :  LMB 
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Figure  C.19  The  Probability  Distributions  of  Cross  Peak  Data  for  FF-1071  Fire 

Pumps  Measured  at  LMB 
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Probability  Probability 


Probability  Density  Function  for  Cross  peak  data 
FF-1062  Fire  Pump  (11/24/89) 

Pickup  Location :  LMB(A) 
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Figure  C.20  The  Probability  Distributions  of  Cross  Peak  Data  for  FF-1062  Fire 

Pumps  Measured  at  LMB(A) 
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Probability  Probability 


Probability  Density  Function  for  Cross  peak  data 
FF-1071  Fire  Pump  (3/10/89) 

Pickup  Location :  LMB(A) 
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Figure  C.21  The  Probability  Distributions  of  Cross  Peak  Data  for  FF-1071  Fire 

Pumps  Measured  at  LMB(A) 
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Probability  Probability 


Probability  Density  Function  for  Cross  peak  data 
FF-1062  Fire  Pump  (11/24/89) 

Pickup  Location :  UPB 
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Figure  C.22  The  Probability  Distributions  of  Cross  Peak  Data  for  FF>1062  Fire 

Pumps  Measured  at  UPB 
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Probability  Probability 


Probability  Density  Function  for  Cross  peak  data 
FF-1071  Fire  Pump  (3/10/89) 

Pickup  Location :  UPB 


(a)  Fire  Pump  #1  (b)  Fire  Pump  #2 
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Figure  C.23  The  Probability  Distributions  of  Cross  Peak  Data  for  FF-1071  Fire 

Pumps  Measured  at  UPB 
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Probability  Probability 


Probability  Density  Function  for  Cross  peak  data 
FF-1062  Fire  Pump  (11/24/89) 

Pickup  Location :  LPB 


(a)  Fire  Pump  #1 
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Figure  C.24  The  Probability  Distributions  of  Cross  Peak  Data  for  FF-1062  Fire 

Pumps  Measured  at  LPB 
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Probability  Probability 


Probability  Density  Function  for  Cross  peak  data 
FF-1071  Fire  Pump  (3/10/89) 

Pickup  Location :  LPB 


(a)  Fire  Pump  #1 
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Figure  C25  The  Probability  Distributions  of  Cross  Peak  Data  for  FF-1071  Fire 

Pumps  Measured  at  LPB 
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APPENDIX  D.  FIGURES  OF  STATISTICAL  ANALYSIS  RESULTS 

This  section  contains  figures  of  the  statistical  analysis  results  with  and  without 
considering  damaged  data  sets. 
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Probability 


Statistical  Analysis  Results 
CG-59  Fire  Pump 
Pickup  Location :  MB(FE) 

Data  analysis  by  using  :  Cross  peak  with  1%  threshold 
Number  of  data  points :  1. 9650^005 
The  sampling  frequency :  1e4004  Hz 


(a)  Linear  Domain  (b)  VdB  Domain 
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Figure  D.l  The  Statistical  Analysis  Results  for  Type  I  Fire  Pumps  Measured  at 
MB(FE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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ProbabiMy 


Statistical  Analysis  Results 
CG-59  Fire  Pump 
Pickup  Location :  MB(CE) 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  2.17304005 
The  sampling  frequency :  1e4004  Hz 


(a)  Linear  Domain 


(b)  VdB  Domain 
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Figure  D.2  The  Statistical  Analysis  Results  for  Type  I  Fire  Pumps  Measured  at 
MB(CE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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ProbabHil 


Statistical  Analysis  Results 
CG-59  Fire  Pump 
Pickup  Location :  PB(CE) 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  3.129e^5 
The  sampling  frequency :  1e^O04  Hz 


(a)  Linear  Domain 
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(b)  VdB  Domain 
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Kurtosis=2.476 


Figure  D.3  The  Statistical  Analysis  Results  for  Type  I  Fire  Pumps  Measured  at 
PB(CE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Rssults 
CG-59  Firs  Pump 
Pickup  Location :  PB(FE) 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  2.614e4005 
The  sampling  frequency :  le^OCH  Hz 


(a)  Linear  Domain 


^=€.002717 
a>::0.002681 
Skewness= 1.253 
Kurtosis^3.959 


(b)  VdB  Donwin 
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Figure  D.4  The  Statistical  Analysis  Results  for  Type  I  Fire  Pumps  Measured  at 
PB(FE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 
CG-59  Rre  Pump 
Pickup  Location :  PB(FE/A) 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  3.105e4005 
The  sampling  frequency :  1e4004  Hz 


(a)  Linear  Domain 


(b)  VdB  Domain 


ii^.003699 
0=0.003307 
Skewness= 1.086 
Kurtosis=3.515 


H  =  106.9 
a  =9.904 
Skewness=‘0.5lS7 
Kurtosis=2.472 


Figure  D.5  The  Statistical  Analysis  Results  for  Type  I  Fire  Pumps  Measured  at 
PB(FE/A)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Probabiiily 


Statistical  Analysis  Results 
DDG  Fire  Pump 
Pickup  Location :  MB(FE) 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  4.733e4004 
The  sampling  frequency :  1e4004  Hz 

(a)  Linear  Domain  (b)  VdB  Domain 


Velocity  (nVs«c)  ^^q-3  VdB  (0VdB-1«-8  nVsec) 
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Figure  D.6  The  Statistical  Analysis  Results  for  Type  n  Fire  Pumps  Measured  at 
MB(FE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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ProtMbiiily 


Statistical  Anaiyais  Raauita 
DDG  Firs  Pump 
Pickup  Location :  MB(CE) 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  4.835e4004 
The  sampling  frequency :  le^OCM  Hz 


(a)  Linear  Domain  (b)  VdB  Domain 
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0=0,0007633 

Skewness=3.243 

Kurtosis=66.02 
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Skewnes5= -0. 0222 
Kurtosis=2.32 


Figure  D.7  The  Stattetical  Analysb  Results  for  Type  II  Fire  Pumps  Measured  at 
MB(CE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Anaiyaia  Raaulta 
DOG  Rra  Pump 
Pickup  Location :  MB(CE/A) 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  pok^ts :  3.089e4004 
The  sampling  frequency :  1e4004  Hz 


(a)  Linear  Domain 


(b)  VdB  Domain 
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Skewness^  1.321 
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H  *  99.98 
a -11.49 
Skewness^O.Ol  939 
Kurtosis- 1.926 


Figure  D.8  The  Statistical  Analysu  Results  for  Type  n  Fire  Pumps  Measured  at 
MB(CE/A)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 
ODG  Fire  Pump 
Pickup  Location :  PB(CE) 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  3.2410^004 
The  sampling  frequency :  104004  Hz 
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Figure  D.9  The  Statistical  Analysis  Results  for  Type  n  Fire  Pumps  Measured  at 
PB(CE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Anaiyaia  Raaulta 
DDG  Fire  Pump 
Pickup  Location :  PB(FE) 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  4.056e4004 
The  sampling  frequency :  1e^<X)4  Hz 


(a)  Linear  Domain 


Velocity  (m/sec) 


(b)  VdB  Domain 


11=0.002166 

0=0.002739 

Skewness=2.008 

Kurtosis=6.446 


H  =  100.8 

o=  10.29 
Skewness= 0. 1 128 
Kurtosis=2.366 


Figure  D.IO  The  Statistical  Analysis  Results  for  Type  II  Fire  Pumps  Measured 
PB(FE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 


Probability 


Statistical  Analysis  Results 
FF-1062  &  FF-1071 
Pickup  Location :  UMB 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  2.1e4005 
The  sampling  frequency :  1e4004  Hz 


(a)  Linear  Domain  (b)  VdB  Domain 


H==0.0006716 

a^.0009258 

Skewness^4.985 

Kurtosis^36.52 


H=‘ 92.12 
a -8.687 
Skewness^O. 02892 
Kunosis-2.953 


Figure  D.ll  The  Statistical  Analysis  Results  for  Type  IH  Fire  Pumps  Measured 
UMB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 
FF-1062  &  FF-1071 
Pickup  Location :  LMB 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  2.176e4005 
The  sampling  frequency :  1e<t<X34  Hz 


(a)  Linear  Domain 


(b)  VdB  Domain 


V 
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0.005  0.01  0.015 

Velocity  (m/sec) 

W  100  150 

VdB  (0VdB*1e-8  nrVsec) 

^=0.0009233 

H- 95.26 

a=0.00l046 

a^S.557 

Skewness-3.35 

Skewness --0.077 39 

Kurtosis=19.59 

Kurtosis-2.697 

Figure  D.12  The  Statistical  Analysis  Results  for  Type  III  Fire  Pumps  Measured  at 
LMB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Statistical  Analysis  Results 
FF-1062  &  FF-1071 
Pickup  Location :  LMB(A) 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  2.Q45e^005 
The  sampling  frequency :  1e+004  Hz 


(a)  Linear  Domain 


Velocity  (nVsec)  ^  ^^-3 


(b)  VdB  Domain 


^l=0.0007025 
a^.000648 
Skewness=  1.881 
Kurtosis=7.538 


93.33 

8.457 

Skewness=  -0. 3486 
Kurtosis=2.6C3 


Figure  D.13  The  Statistical  Analysis  Results  for  Type  ID  Fire  Pumps  Measured 
LMB(A)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 


at 
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Statistical  Analysis  Results 
FF-1062  &  FF-1071 
Pickup  Location :  UPB 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  1.367e4005 
The  sampling  frequency :  le>004  Hz 


Qi  II  U  11  H  U— I - -  I 

0  0.02  0.04 

Velocity  (m/sec) 


^=0.004141 

<7=0.004442 

Skewness=1.57 

Kunosis=5.446 


(b)  VdB  Domain 


VdB  (OVdBsle-S  rrVsec) 

11  =  106.1 
<7=11.84 
Skewness= -0.4083 
Kurtosis=2.207 


Figure  D.14  The  Statistical  Analysis  Results  for  Type  III  Fire  Pumps  Measured  at 
UPB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Probability 


Statistical  Analysis  Results 
FF-1062  &  FF-1071 
Pickup  Location :  LPB 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points ;  1.99e>005 
The  sampling  frequency :  Hz 


(a)  Linear  Domain 


(b)  VdB  Domain 


0  0.01  0.02  0.03 

Velocity  (m/sec) 


_ I  ■  1  a  ■■  ■  II  «  ■  r-i _ I 

50  100  150 

VdB  (0VdBs1e-8  m/sec) 


H=0.00316 
a^.003246 
Skewness=l.  779 
Kurtosis-6.775 


105.1 
a  =10.2 

Skewness=  -0. 3726 
Kurtosis=2.43 


Figure  D.15  The  Statistical  Analysis  Results  for  Type  in  Fire  Pumps  Measured  at 
LPB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Probability 


Statistical  Analysis  Results 
CG-59  Fire  Pump 

Pickup  Location  :  PB(CE)  w/o  89#6,90#3,90#6 

Data  e  alysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  2.704e4005 
The  sampling  frequency :  1e+004  Hz 


(a)  Linear  Domain 


Velocity  (nVsec) 

^=0.002319 
<7=0.001988 
Skewness^  1.046 
Kurtosis=3.304 


(b)  VdB  Domain 


VdB  (0VdB>1e^  nrVsec) 

11  =  103.3 
0  =  9.245 
Sk£yvness='0.5101 
Kurtosis=2.506 


Figure  D.16  The  Statistical  Analysis  Results  (Without  Damaged)  for  Type  I  Fire 
Pumps  Measured  at  PB(CE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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ProbabHily 


Statistical  Analysis  Results 
CG-59  Fire  Pump 

Pickup  Location  :  PB(FE)  w/o  90#3 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  2.489e4005 
The  sampling  frequency ;  1e>f004  Hz 

(a)  Linear  Domain  (b)  VdB  Domain 


/i=0.002525 
a^.002472 
Skewness^  L 169 
Kurtosis^3.409 


yL^lOS 
a -10.3 

Skewness9-0.3165 

Kurtosis^2.228 


Figure  D.17  The  Statistical  Analysis  Results  (Without  Damaged)  for  Type  I  Fire 
Pumps  Measured  at  PB(F£)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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ProbabHily 


Statistical  Analysis  Results 
DDG-994  Fire  Pump 
Pickup  Location :  MB(CE)  w/o  #1 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  3.238e-i<X)4 
The  sampling  frequency :  1e4004  Hz 


(a)  Linear  Domain 


Velocity  (m/sec) 


(b)  VdB  Domain 


H=0.0006192 

(P=0.0008105 

Skewness=2.l46 

Kunosis-6.762 


=  90.52 
0=9.194 
Skewness=0. 4951 
Kurtosis=2.566 


Figure  D.18  The  Statistical  Analysis  Results  (Without  Damaged)  for  Type  II  Fire 
Pumps  Measured  at  MB(CE)  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Probabifily 


Statistical  Analysis  Results 
FF.1062  &  FF-1071 
Pickup  Location  :  UMB  w/o  62#1 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  2.05e4005 
The  sampling  frequency :  1e4004  Hz 


(a)  Linear  Domain 


Velocity  (m/sec)  ^^^-3 


(b)  VdB  Domain 


60  80  100  120 
VdB  (OVOBsle-S  m/sec) 


H=0.0005973 
a^.0006017 
Skewness^2.488 
Kurtosis^  11.46 


fi  =  91.81 
0^8.338 
Skewness^-0. 1669 
Kurtosis^2.595 


Figure  D.19  The  Statistical  Analysis  Results  (Without  Damaged)  for  Type  III  Fire 
Pumps  Measured  at  UMB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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ProbabiMy 


Statistical  Analysis  Results 
FF-1062  &  FF-1071 

Pickup  Location  :  LMB  w/o  62#1,71#2 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  1.873e4005 
The  sampling  frequency :  1ef004  Hz 


(a)  Linear  Domain 


(b)  VdB  Domain 


0  0.005  0.01 

Velocity  (m/sec) 


H=0.0007757 
a^.0007364 
Skewness— 1.98 
KurtosL^  7.43 


60  80  100  120 
VdB  (0VdB-1e^  m/$ec) 


11  =  94.27 


c=8.15 

Skewness='0.181 

Kurtosis=2.555 


Figure  D.20  The  Statistical  Analysis  Results  (Without  Damaged)  for  Type  III  Fire 
Pumps  Measured  at  LMB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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Probability 


Statistical  Analysis  Results 
FF-1062  &  FF-1071 

Pickup  Location :  UPB  w/o  FF-1062  #1 

Data  analysis  by  using :  Cross  peak  with  1%  threshold 
Number  of  data  points :  1 .2450+005 
The  sampling  frequency :  1e+004  Hz 


(a)  Linear  Domain 


(b)  VdB  Domain 


0.01  0.02 
Velocity  (rrVsoc) 


H=0.003603 
<j=0.003906 
Skewness- L7  11 
Kurtosis^6.361 


Q'  ■...■■■■..Ml  -1 
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VdB  (0V(ffi>1e-8  m/sec) 


105.1 


Skewness^  -0.3666 
Kurtosis=2.176 


Figure  D.21  The  Statistical  Analysis  Results  (Without  Damaged)  for  Type  ni  Fire 
Pumps  Measured  at  UPB  (a)  Linear  Velocity  Domain  (b)  VdB  Domain 
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APPENDIX  E.  FIGURES  OF  1/1  OCTAVE  BAND  LEVELS 

This  section  contains  all  the  figures  of  1/1  octave  band  levels. 
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VdB  VdB  VdB 


1/1  Octave  Band  Analyaia 
CG-59  Fire  Pump  (10/30/89) 
Pickup  Location :  MB(FE) 

0  VdBsle*8  m/sec 


(a)  Pump  #1 

10-5k  Hz  B.B.  L*v«l>  97^1  dB . 


(b)  Punp  #2 

10-5k  Hz  B.B.  Lavakt  95.52  dB . 
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(c)  Pump  #3 

20  1 0*Sk  Hz  B.B.  Lavals  93.99  dB  . 
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(a)  Pump  #5 

^20 [  10-5k  Hz  B.B.  Lavab  98.70  dB  . 


(f)  Pump  #6 

,2o[  10-5k  Hz  B.B.  Lavak  96.76  dB . 


B  100 


10*  10*  10* 
Fraquancy  (Hz) 


10’  10*  10*  10* 
Fraquancy  (Hz) 


Figure  E.1  The  1/1  Octave  Band  Analysis  Results  for  CG-59(10/30/89)  Fire  Pumps 

Measured  at  MB(FE) 
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VdB  VdB  VdB 


1/1  Octave  Band  Analysis 
CG-59  Fire  Pump  (2/26/90) 
Pickup  Location :  MB(FE) 

0  VdBs1e*8  m/sec 


(a)  Pump  #1 

,2of  10-5k  Hz  B.B.  L«v«b  96.17  dB 


(b)  Pump  #2 

,2o[  I^Sk  B.B.  Lavak  95.65  dB . 
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(e)  Pump  #5 


120  10-5k Hz  B.B.  Levels 97.97 dB. 
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(f)  Pump  #6 

120 [  10-5k  Hz  B.B.  Levels  96.97  dB . 
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Figure  E.2  The  1/1  Octave  Band  Analysis  Results  for  CG>59(2/26/90)  Fire  Pumps 

Measured  at  MB(FE) 


142 


VdB  VdB  VdB 


1/1  Octave  Band  Analyais 
CG-59  Fire  Pump  (10/30/89) 
Pickup  Location :  MB(CE) 

0  VdBsle-8  m/sec 


(a)  Pump  «1 

^20  10-5k  Hz  B.B.  L«v«t>  99.18  dB^ 


(b)  Pump  #2 

10-5k  Hz  B.B.  Laval-  9929  dB . 
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(c)  Pump  #3 
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(e)  Pump  #5 

120  [  10-5k  Hz  B.B.  Leveb  98.90  dB . 


(()  Pump  #6 

120  [  10-5k  Hz  B.B.  Level-  9824  dB . 
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Figure  E  J  The  1/1  Octave  Band  Analysis  Results  for  CG-S9(1 0/30/89)  Fire  Pumps 

Measured  at  MB(CE) 
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VdB  VdB  VdB 


1/1  Octave  Band  Analysis 
CG-59  Fire  Pump  (2/26/90) 
Pickup  Location :  MB(CE) 

0  VdBs1e-8  m/sec 
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Figure  E.4  The  1/1  Octave  Band  Analysis  Results  for  CG-59(2/26/90)  Fire  Pumps 

Measured  at  MB(CE) 
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VdB  VdB  VdB 


1/1  Octave  Band  Analysis 
CG-59  Fire  Pump  (10/30/89) 
Pickup  Location :  PB(CE) 

0  VdB=1e-8  m/sec 


(a)  Pump  «1 

,2of  10-5k  Hz  B.B.  L«ve»=107.77  dB . 


(b)  Pump  «2 

^20  10-5k  Hz  B.B.  Lavabl 06.64 dB. 
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(c)  Pump  #3 

,2o[  10-5k  Hz  B.B.  Level=109.41  dB  . 


(d)  Pump  #4 
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(e)  Pump  #5 

^20  10-5k  Hz  B.B.  Levebi  07.83  dB . 
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Figure  £.5  The  1/1  Octave  Band  Analysis  Results  for  CG-59(10/30/89)  Fire  Pumps 

Measured  at  PB(CE) 
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VdB  VdB  VdB 


1/1  Octave  Band  Analysis 
CG-59  Fire  Pump  (2/26/90) 
Pickup  Location :  PB(CE) 

0  VdBsile-S  m/sec 


(a)  Pump  #1 

t0-5k  Hz  B.6.  UvAkloa.ei  dB . 
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.,2o|  10-5k  Hz  B.B.  LeveM 07.06  dB. 
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(d)  Pump  #4 
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Figure  £.6  The  1/1  Octave  Band  Analysis  Results  for  CG-59(2/26/90)  Fire  Pumps 

Measured  at  PB(CE) 
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VdB  VdB  VdB 


1/1  Octave  Band  Analysis 
CG-59  Fire  Pump  (10/30/89) 
Pickup  Location :  PB(FE) 

0  VdB*1e-8  m/sec 


(a)  Pump  #1 

^20  1 0-5k  H2  B.B.  L«v«is1 09^8  dB . 


(b)  Pump  *2 

10-Sk  Hz  B.B.  Lavak103.30  dB . 
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(c)  Pump  #3 

,2o[  10-5k  Hz  B.B.  Leveb107.15  dB . 


(d)  Pump  «4 

^2o[  10-5kHzB.B.  Leveb109.16dB. 
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(e)  Pump  #5 

^20  10-5k  Hz  B.B.  Leveb109.07  dB . 


(f)  Pump  #6 

.,2o[  10-5k  Hz  B.B.  Levebi  09.58  dB . 
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Figure  E.7  The  1/1  Octave  Band  Analysis  Results  for  CG-59(1 0/30/89)  Fire  Pumps 

Measured  at  PB(FE) 
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VdB  VdB  VdB 


1/1  Octave  Band  Analysis 
CG-59  Fire  Pump  (2/26/90) 
Pickup  Location :  PB(FE) 

0  VdBs1e*8  m/sec 
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Figure  £.8  The  1/1  Octave  Band  Analysis  Results  for  CG-59(2/26/90)  Fire  Pumps 

Measured  at  PB(FE) 
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VdB  VdB  VdB 


1/1  Octave  Band  Analyaia 
CG-59  Fire  Pump  (10/30/89) 
Pickup  Location :  PB(FE/A) 

0  VdBs1e-8  m/sec 


(a)  Pump  «1 

^20 [  10-5k  Hz  B.B.  L«v«lx112.17  dB . 


(b)  Pump  «2 
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^20 1  10-5k  Hz  B.B.  Level=1 1 1 .95  dB  . 


(0  Pump  #6 

^2o[  l0-5k Hz  B.B.  LeveM 09.68 dB. 


§  100 


10*  10*  10* 
Frequency  (Hz) 


10’  10*  10*  10* 
Frequency  (Hz) 


Figure  E.9  The  1/1  Octave  Band  Analysis  Results  for  CG*59(10/30/89)  Fire  Pumps 

Measured  at  PB(FE/A) 
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SPA  SPA  SPA 


1/1  Octave  Band  Analyaia 
CG-59  Fire  Pump  (2/26/90) 
Pickup  Location :  PB(FE/A) 

0  VdBs1e-8  m/sec 


(a)  Pump  *1 

10-5k  Hz  B.B.  L«v«Ib108.93  dS 
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Figure  E.10  The  I /I  Octave  Band  Analysis  Results  for  CG-59(2/26/90)  Fire  Pumps 

Measured  at  PB(FE/A) 
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VdB  VdB 


1/1  Octavo  Band  Analyaia 
DDG-994  Fire  Pump 
Pickup  Location :  MB(FE) 


0  Vd6»1e-8  m/sec 


(a)  Pump  #1 


Frequancy  (Hz) 


(b)  Pump  *2 


(c)  Pump  #3 


Figure  £.11  The  1/1  Octave  Band  Analysis  Results  for  DDG>994  Fire  Pumps 

Measured  at  MB(FE) 
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1/1  Octave  Band  Analysis 
DDG-994  Fire  Pump 
Pickup  Location :  MB(CE) 
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Figure  £.12  The  1/1  Octave  Band  Analysis  Results  for  DDG-994  Fire  Pumps 

Measured  at  MB(CE) 
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1/1  Octave  Band  Analyaia 
DDG-994  Fire  Pump 
Pickup  Location :  MB(CE/A) 
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(c)  Pump  #3 


Figure  E.13  The  1/1  Octave  Band  Analysis  Results  for  DDG-994  Fire  Pumps 

Measured  at  MB(CE/A) 
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1/1  Octave  Band  Analysis 
DDQ-994  Fire  Pump 
Pickup  Location :  PB(CE) 
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Figure  £.14  The  1/1  Octave  Band  Analysis  Results  for  DDG-994  Fire  Pumps 

Measured  at  PB(CE) 
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1/1  Octave  Band  Analyaia 
DDG-994  Fire  Pump 
Pickup  Location :  PB(FE) 

0  VdBs1e-8  m/sec 

(a)  Pump  #1  (b)  Pump  #2 


Frequency  (Hz)  Frequency  (Hz) 


(c)  Pump  #3 


Figure  E.15  The  1/1  Octave  Band  Analysis  Results  for  DDG-994  Fire  Pumps 

Measured  at  PB(FE) 
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1/1  Octave  Band  Analysis 
FF-1062  Fire  Pump  (11/24/89) 
Pickup  Location :  UMB 
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Figure  E.16  The  1/1  Octave  Band  Analysis  Results  for  FF-1062  (11/24/89)  Fire  Pumps 

Measured  at  UMB 
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1/1  Octave  Band  Analysis 
FF-1071  Fire  Pump  (3/10/89) 
Pickup  Location :  UMB 
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Figure  E.17  The  1/1  Octave  Band  Analysis  Results  for  FF-1071(3/10/89)  Fire  Pumps 

Measured  at  UMB 
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1/1  Octave  Band  Analysis 
FF-1062  Fire  Pump  (11/24/89) 
Pickup  Location :  LMB 
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Figure  £.18  The  1/1  Octave  Band  Analysis  Results  for  FF-1062  (1 1/24/89)  Firt  Pumps 

Measured  at  LMB 
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1/1  Octave  Band  Analysis 
FF-1071  Fire  Pump  (3/10/89) 
Pickup  Location :  LMB 
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Figure  E.19  The  1/1  Octave  Band  Analysis  Results  for  FF>1071(3/10/89)  Fire  Pumps 

Measured  at  LMB 
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1/1  Octave  Band  Analysis 
FF-1062  Fire  Pump  (11/24/89) 
Pickup  Location :  LMB(A) 
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Figure  E.20  The  1/1  Octave  Band  Analysis  Results  for  FF-1062  (11/24/89)  Fire  Pumps 

Measured  at  LMB(A) 
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1/1  Octave  Band  Analysis 
FF-1071  Fire  Pump  (3/1Q/89) 
Pickup  Location :  LMB(A) 
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Figure  E.21  The  1/1  Octave  Band  Analysis  Results  for  FF-1071(3/10/89)  Fire  Pumps 

Measured  at  LMB(A) 
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1/1  Octave  Band  Analysis 
FF-1062  Fire  Pump  (11/24/89) 
Pickup  Location :  UPB 
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Figure  £.22  The  1/1  Octave  Band  Analysis  Results  for  FF-1062  (11/24/89)  Fire  Pumps 

Measured  at  UPB 
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Figure  E.23  The  1/1  Octave  Band  Analysis  Results  for  FF- 107 1(3/1 0/89)  Fire  Pumps 

Measured  at  UPB 
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1/1  Octave  Band  Analyais 
FF-1062  Fire  Pump  (11/24/89) 
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Figure  E.24  The  1/1  Octave  Band  Analysis  Results  for  FF-1062  (11/24/89)  Fire  Pumps 

Measured  at  LPB 
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1/1  Octave  Band  Analyaia 
FF-1071  Fire  Pump  (3/1Q/89) 
Pickup  Location :  LPB 


0  VdB«le-8  n/sec 
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(c)  Pump  #3  (d)  Pump  #4 
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Figure  E.25  The  1/1  Octave  Band  Analysis  Results  for  FF-1071(3/10/89)  Fire  Pumps 

Measured  at  LPB 
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APPENDIX  F.  FIGURES  OF  1/1  OCTAVE  BAND  LEVELS  USING 
NAVSSES  FREQUENCY  DOMAIN  DATABASE 


This  section  contains  all  the  figures  of  1/1  octave  band  levels  using  NAVSSES 
frequency  domain  database. 
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Figure  F.l  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (05/12/93)  Fire  Pump#l 
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Figure  F.2  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (05/12/93)  Fire  Punip#2 
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Figure  F.3  The  1/1  Octave  Band  Analysis  Results  for  CG*S9  (05/12/93)  Fire  Punip#3 
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Figure  F.4  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (05/12/93)  Fire  Punip#4 
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Figure  F>5  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (05/12/93)  Fire  Punip#5 


171 


1/1  Octave  Band  Analysis 
CG-59  Fire  Pump#6 
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Figure  F.6  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (05/12/93)  Fire  Pump^ 
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1/1  Octave  Band  Analysis 
CG-59  Fire  Pump#1 
(02/23/92) 
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Figure  F.7  The  I/l  Octave  Band  Analysis  Results  for  CG-59  (02/23/92)  Fire  Pump#! 


173 


VdB  VdB 


1/1  Octave  Band  Analysis 

- 

120 

CG*59  Fire  Pump#2 

(02/23/92) 

O  100 

5 

80 

0  VdB=1e-8  m/sec 

60 

10 

(a)  MB(FE) 

B.B.  L«v«l=106.41  dB 


10*  10*  10* 
Fraquancy  (Hz) 


(b)  MB(CE) 

B.B.  Lavel=107.69  dB 


(c)  PB(CE) 

B.B.  Levels lie.lOdB 


m  100 

5 


10*  10*  10* 
Frequency  (Hz) 


10*  10*  10* 
Frequency  (Hz) 


(d)  PB(FE) 

B.B.  Level=116.32dB 


(e)  PB(FE/A) 

B.B.  Level=119.39dB 


CD  100 

5 


10‘  10* 
Frequency  (Hz) 


10*  10*  10* 
Frequency  (Hz) 


Figure  F.8  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (02/23/92)  Fire  Punip^2 
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Figure  F.9  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (02/23/92)  Fire  Punip#3 
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1/1  Octave  Band  Analysis 
CG-59  Fire  Pump#4 
(02/23/92) 
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Figure  F.IO  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (02/23/92)  Fire  Punip#4 
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1/1  Octave  Band  Analysis 
CG-59  Fire  Pump#5 
(02/23/92) 
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Figure  F.ll  The  1/1  Octave  Band  Analysis  Results  for  CG>59  (02/23/92)  Fire  Pump^S 
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Figure  F.12  The  1/1  Octave  Band  Analysis  Results  for  CG-S9  (10/21/90)  Fire  Pump^l 
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Figure  F.13  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (10/21/90)  Fire  Pump#2 
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Figure  F.14  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (10/21/90)  Fire  Punip#3 
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Figure  F.15  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (10/21/90)  Fire  Puniip#4 
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Figure  F.16  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (10/21/90)  Fire  Pump#S 
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Figure  F.17  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (06/07/90)  Fire  Pump#! 
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1/1  Octave  Band  Analysis 
CG-59  Fire  Pump#2 
(06/07/90) 
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Figure  F.18  The  1/1  Octave  Band  Analysis  Results  Tor  CG-59  (06/07/90)  Fire  Punip#2 
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Figure  F.19  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (06/07/90)  Fire  Punip#3 
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Figure  F.20  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (06/07/90)  Fire  Puinp#4 
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Figure  F.21  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (02/26/90)  Fire  Pump^l 
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Figure  F.22  The  1/1  Octave  Band  Analysis  Results  for  CG>59  (10/30/89)  Fire  Pump#! 
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Figure  F.23  The  1/1  Octave  Band  Analysis  Results  for  CG-S9  (10/30/89)  Fire  Punip#2 
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Figure  F.24  The  1/1  Octave  Band  Analysis  Results  for  CG>S9  (10/30/89)  Fire  Punip#3 
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Figure  F.25  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (10/30/89)  Fire  Punip#4 
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Figure  F.26  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (10/30/89)  Fire  Punip#5 


192 


1/1  Octave  Band  Analysis 
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Figure  F.27  The  1/1  Octave  Band  Analysis  Results  for  CG-59  (10/30/89)  Fire  Punip#6 
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